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Abstract  of  Dissertation  Presented  to  the  Graduate  School 
of  the  University  of  Florida  in  Partial  Fulfillment  of  the 
Requirements  for  the  Degree  of  Doctor  of  Philosophy 

ON  THE  DISTRIBUTION  OF  THE  LENGTH  OF  A 
SPHERICAL  RANDOM  VECTOR 

BY 

Everton  De  Courcey  Rowe 
August  1988 

Chairman:  John  G.  Saw 

Major  Department:  Statistics 

This  dissertation  is  concerned  with  determining  the 
probability  density  function  of  the  length  of  a spherical 
random  vector. 

The  case  of  the  sum  of  independent  identically 
distributed  random  vectors,  each  of  which  has  a uniform 
distribution  on  the  unit  m-sphere  is  first  investigated. 
Results  which  are  similar  to  those  for  the  three 
dimensional  situation  are  established.  The  probability 
density  function  of  the  length  of  the  sum  is  given  in  terms 
of  the  probability  density  function  of  the  absolute  value 
of  a single  element  of  that  vector. 

The  next  case  is  to  investigate  the  general  family  of 
spherical  densities.  We  do  so  by  two  methods.  First  we  use 
the  characteristic  function  to  develop  the  main  result  for 
an  odd  number  of  dimensions.  Then  we  derive  an  integral 


-vi  - 


equation  which  is  useful  tor  all  values  of  m (m>2)  , the 
number  of  dimensions.  We  use  this  integral  equation  to 
reproduce  some  known  results.  It  is  made  clear  that  the 
results  established  for  the  uniform  distribution  on  the 
unit  sphere  is  a special  case  of  this  more  general  result. 
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CHAPTER  I 
INTRODUCTION 


1 . 1 Directional  Data 


quite  often  scientific  data,  are  recorded  as  unit 


vectors . 

As 

such  , 

these  vectors  can 

be 

considered  as 

d i rect i ons 

i n 

space 

and  are  referred 

to 

as  directional 

data.  In  two  dimensions  these  vectors  correspond  to  points 
on  the  unit  circle.  In  three  dimensions  these  points  are 
located  on  the  surface  of  a sphere.  For  higher  dimensions 
we  speak  of  these  points  as  lying  on  the  surface  of  a 
hypersphere . 

The  study  of  directional  data  was  first  focused  on  the 
two  dimensional  and  three  dimensional  situation.  However, 
one  can  derive  results  for  a general  number  of  dimensions 
using  the  same  techniques  which  are  required  for  the  study 
of  two  or  three  dimensions;  this  is  easily  seen  in  the  case 
of  a random  vector  which  has  a uniform  distribution  on  the 
sphere . 

Results  pertaining  to  higher  than  three  dimensions  are 
of  practical  importance,  since  data  may  be  collected  on 
several  variables.  A vector  corresponding  to  information 
on  each  of  these  variables  may  be  scaled  so  that  the  length 
of  the  vector  is  unity. 
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The  usual  methods  for  analyzing  multivariate  data  are 
based  on  the  multivariate  Gaussian  distribution,  but  these 
methods  are  not  applicable  here  and  so  an  investigation  of 
distributions  on  the  sphere  is  of  interest  to 
st at  i st  i c i ans  . 

Mardia  (1972)  is  a good  introduction  to  the  literature 
(prior  to  1972)  on  directional  data.  Watson  (1982)  deals  in 
detail  with  a random  vector  which  has  a uniform 
distribution  on  the  unit  sphere.  This  book  provides  an 
indication  of  the  direction  in  which  research  in  this  area 
is  going  and  it  references  many  of  the  recent  contributions 
to  the  subject. 

Directional  data  have  been  useful  in  diverse  areas. 
Amongst  these  are  sampling,  astronomy,  medicine, 
psychology,  physics,  biology,  geography  and  meteorology. 
Some  of  the  main  models  used  to  analyze  directional  data 
are  the  Bingham  distribution,  the  uniform  distribution  on 
the  unit  sphere,  the  Langevin  distribution,  Saw’s 
distribution  and  the  von  Mises  distribution  on  the  circle. 

In  this  chapter  we  study  the  uniform  distribution  on 
the  unit  m-sphere . We  consider  an  arbitrary  number  of 
dimensions.  In  Section  1.2  we  introduce  the  notation  which 
we  use  throughout  this  chapter.  Section  1.3  contains  some 
distributional  results  for  a random  vector  which  has  a 
uniform  distribution  on  the  unit  m-sphere.  Section  1.4  is 
a summary  of  a paper  by  Saw  (1986) . 
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1 . 2 Notat  i on 


In  what  follows  Rm  denotes  m-d imens i onal  Euclidean 
space . 

i 

With  || x | = (x'x)  2 let  Sm={x  e Rm : |x||  = 1}  be  the  unit  m- 
sphere  and  Bm  = {x.  £ Rm:|x||  < 1}  be  the  unit  m-bal  1 . 

It  W is  a random  variable  with  probability  density  f (w) 
given  by 


f (w) 


(1-w2)  2 

0, 


-1<w<1 , 


e 1 sewhere , 


we  will  write  W~WI/.  Throughout,  we  let  v = 
B(p,q)  = is  the  usual  beta 

r(p+q) 

P(X<t)=P(Y<t)  V t e R1  we  write  X~Y. 


variable  having  a probability  density  given 


im-1 . Here 
function.  If 

For  a random 
by 


SO) 


x^Cl-x)6"1 

B (a , b) 


0 


0<x<  1 


e 1 sewhere , 


we  write  X~B(a,b). 

Wilks  (1962)  gives  the  gamma  distribution  as  a one 
parameter  family;  we  adopt  the  same  convention  here.  Thus 
if  Y has  a gamma  density,  then  the  probability  density 
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function  is  given  by 


g(y) 


u—\  — y 
y e 

I>) 


0<y <oo , 


0 elsewhere, 


and  we 

write  Y~ 

GO)  . 

As  usual , Xm  denotes  a 

ch  i 

square 

random 

variable 

with 

m degrees  of  freedom  and 

Xm 

i s the 

positive  square  root  of  a \m  random  variable.  If 

* (Yu  > y*)  = 


..  v 1-1 

f(^i)  r(r2)...ro*)yi 


(l-yj-  . . . -yt) 


9 


where  0<yj<l  , i=l,2,...,k  and  0<yi+y2+ . . . +yi.<  1 . We  say 

yi  »y2  > • • • >y*  have  a Dirichlet  distribution,  denoted  by 
DCi/j,^,  . . . , vk  : i/i+1)  . 

If  a random  vector  Y has  a multivariate  normal 
distribution  with  mean  vector  ji  and  variance-covariance 
matrix  , we  write  this  as  Y~N(^£,  JZ)  . 

We  define  the  function  Al/(t)  by  the  infinite  series 


Aj/  (t ) 


=£ 


k=0 


/i1-2n  k r(r  + l) 

' k ! T (r+k+1 ) 


V fit2!1 1 

t=o  2 k!  m(m+2) . . . (m+2k-2)' 


The  function  A^t)  is  related  to  the  Bessel  function 


Jl/(t)  and  the  modified  Bessel  function  Ij/(t)  by 
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Ir  (t) 


and 


h-u  (t ) 

r (t'+i ) ’ 


Jr  (t ) 


a*y 


A!/ (it) 

r(v+i) ' 


The  coefficient  oT  a n,  in  the  power  series  expansion  oT 
the  function  ( 1 -2at+a2)  ~v , is  Cn(t)  and  is  called 
Gegenbauer  polynomial  of  degree  n.  The  fact  that  Co(t)=l 
will  be  useful.  Finally,  the  confluent  hypergeometric 
function  M(a,b:z)  is  defined  by  the  infinite  series 


M (a , 


oo 

b:z)  = Yj 


*=o 


z*(a)fc 

k!(b)fc’ 


where  (a)j.  in  Pochammer’s  symbol  and  is  a(a+l)  . . . (a+k-1) 
and  (a)  0 = 1 . 


1*3  The  Uniform  Distribution  on  the  Unit  Sphere 


For  -l<t<l  and  ( e Sm , let  Cm((,t)={u  e Sm  :(fu>t}  and 
A be  the  smallest  <r-algebra  generated  by  }Cm((,t):  (eSm, 
t £ [-1,1]};  thus, 

A = cr  {Cm(C_,t):C_  £ Sm , t £ [-1,1]}. 

Furthermore  let  Um  be  the  measure  given  by 

ii  rn  (r  fU-  Surface  area  of  Cm(C,t) 
m\  m(s>  )}  Surface  area  of  Sm 
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Transf ormation  by  multiplying  by  an  orthogonal  matrix 
is  a rotation  of  the  original  figure  and  thus  areas  and 
angles  in  the  original  figure  are  preserved  under  such  a 
transformation.  It  follows  that  for  A e A and  for  P,  an 
orthogonal  matrix,  we  have  Um (A)  = Um(PA),  and  consequently 
Um (Cm (C  > t) ) does  not  depend  on  £.  The  triplet  (Sm  , A , Um) 
is  a probability  space.  We  can  construct  any  number  of 
densities  on  the  unit  sphere  Sm  with  respect  to  the  measure 
Um . However  in  this  dissertation  we  restrict  our  attention 
mainly  to  the  density 

f(u)  =1  U € Srn  , 

we  say  U has  a uniform  distribution  on  the  unit  m-sphere 
and  write  U ~U(Sm). 

We  can  now  give  a characterization  of  a random  vector  U 
having  a uniform  distribution  on  the  unit  sphere. 

Theorem  1 . Let  Z be  an  mxl  component  vector  of  independent 
identically  distributed  Gaussian  random  variables  with 
§ Zj=  0 and  § Zj=  l.  Then 

(i)  |firu(Sm)’ 

(ii)  pry,  and  ||Z||  are  independent,  and 

IzHI 

II 2 2 

A | ~Xm • 


(iii) 
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Proof : Make  the  transformation  R=|Z||  and  U = r^rr,  ■ The 

m 

Jacobian  of  “this  ‘transformation  is  R"1"*1 . Thus 


r 


m — 1 


f(r,u)  _ r 

nif™-’) 


_Ir2 
. 2 


ram) 


r(» 

Kr(i))r 


0<r<oo , |U||  = 1 


which  is  the  result  we  set  out  to  prove. 


This  theorem  is  very  useful . We  can  obtain  results 

about  the  components  of  a random  vector  U having  a uniform 

distribution  of  the  unit  sphere  by  examining  the 

corresponding  components  of  a vector  where  Z~N(0,I). 

U±\ 

Our  familiarity  with  Z greatly  simplifies  the  task. 


Result  1 . For  every  fixed  mxm  orthogonal  matrix  P,  we  have 
PU~U (Sm)  . 


Proof : 
f ol lows 


Theorem 
that  PU~ 


tell  us  that 
. Making  the 


LJ~|§t,  where  Z~N(0,I). 
— 

transformation  Y = PZ 


It 

and 


noting  that 
Y~N (0,1) . 


PZ 

|Z|| 


PZ 


J(Z'p'PZ) 


Y 

II Y II 


we  have  PU  '■ 


whf 


Result  2.  For  every  fixed  unit  vector  a,  a,U~WI/ . 


Proof : Without  loss  of  generality  we  can  consider  a'  as 
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the  first  row  of  the  matrix  P in  Result  1 

/ 2 

ii2  ' ~ “V2’2 


It  is  now 


i Y Y2 

obvious  that  a U~  |yj . We  know  thatjp^i  ~ B(|,|(m-i)).  The 


result  is  immediate  on  making  the  square  root 
transformation  and  on  noting  that 


1X8  II  y|  ’ 

We  now  give  a more  general  version  of  Result  2;  but  we 
first  give  as  lemmas,  two  results  from  Wilks. 


Lemma  1 . 


Suppose  X!,X2,  . . . independent  rand 


om 


variables  having  gamma  distributions  G(iq)  , . . . ,G(i/j.+1)  . 
Make  the  transformation, 


Y.-  = 


X, 


X^+Xj-j- . . . + X 


jfc+1 


i=l ,2, . . . ,k, 


then  Yj  , Y2,...,Yt  has  the  k-variate  Dirichlet  distribution 
D(i/1,i/2,  . . . ,vk:  vk+l)  . 

Lemma  2 . If  (Xl5X2,  . . . ,Xfc)  is  a random  vector  having  the  k- 
variate  Dirichlet  distribution  then  the  marginal 

distribution  of  (Xx  ,X2 , . . . ,Xt  ) , k^k  is  the  k-variate 
Dirichlet  distribution  D (iq  , . . . , vk^ : vk^y+  . . .+rjfe^_1)  . 

Result  3.  Partition  U/  as  (U,1,U,2)  where  U,i=(U1,U2,  . . . ,Uj.) 
and  = (Ut+1 , . . . ,Um)  then 


(i) 


Hi 

lu.l 


u(sk) 
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(ii)  || y 1 1|  and  || JJ x ||  are  independent, 


(iii)  \\Uif  ~B(fk,  i(m-k))  and 


(iv)  The  probability  density  of  Ux , say  g^i)  , is  given  by 


g(u.i) 


r(» 

(r(l))kr(I(m-k)) 


(i  - n/ui) 


ik-i 


0 < u/ni  < 1, 


Proof : Equation  (i)  is  immediate  on  noting  that 

|Z|  _ z , 

|Ui|  ||Z||  • IZJ  - IZJ  ’ 

and  by  Theorem  1 we  obtain  the  result  that  „ ~~1  „ has  a 

ll^il 

uniform  distribution  on  the  unit  k-sphere. 


Now 


Hi 


z z 

ll ~l  ||  and  by  Theorem  1 , we  have  that  — 1 


m 


z 

is  independent  of  |21|.  It  follows  that  n"1,  is  independent 


of  ||Z|| , and  consequently 


Zi  . 


m 


is  independent  of 


I u 
m 


wh  i ch 


is  a restatement  of  (ii) . 


Equat i on ( i i i ) is  immediate  on  noting  that 
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|Ui|2' 


£z  ? 

i — l 


Ez  ? + £zf 

i 


fc  + 1 


where  Z ~ N (0,  I).  It  is  clear  that  Z?  i = l , 2 , . . . , m are 

independent  %i  random  variables.  In  the  notation  of  Wilks 


lZ?  ~ G(l),  and  on  using  Lemma  1 it  is  clear  that 


_Zj 

llzll  ’ 


i = 1 , . . . ,m-l  have  a D(1,1,...,L:1)  distribution . Lemma  2 

Z • 

tells  us  that  the  random  variables  j—4  i=l,...,k  have  a 

|Z| 

Dirichlet  distribution  D (1 , 1 , . . . , 1;  On  making  the 

symmetric  square  root  transformation  and  on  multiplying  by 


the  Jacobian  of  this  transformation  we  see  that 
density  given  in  (iv) . 


|Zi| 


has  the 


Result  4.  Let  U be  a vector  having  a uniform  distribution 
on  the  unit  sphere,  then 

(i)  g U =0  and 

(ii)  g u u'  = i I. 

Proof : For  j=  1,2,. . . ,m  it  is  clear  that 


°=  6 Z,-  = g -^jljp.IZfl  = g -j|j-  g \U  . 


It  follows  that 
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e u>-  = s IZI 


J _ 


= 0 


Furthermore 

««  = $ z<z/  = 8 I • t • R 


£ £ r2 

“ S R R ' S H ’ 


i=J 


where  < — 


0 


i#j 


5,--  „ 

Hence,  § U,Uj  = - — — ^ and  8 R = m,  and  the  proof  is  complete 

8 R 


From  the  multivariate  central  limit  theorem,  Rao  (1965) 
section  2c. 5,  it  is  clear  that  it  U,  i = l , . . . , n are 
independent  and  identically  distributed  random  vectors  with 
U1~U(Stn)  , then 

•\|(m/n)  £u,  4 N (0 , I) 

i=i 

and 

(m/n)  ^ X2m 

Result  5.  The  characteristic  Function  <^>(t),  and  the  moment 
generating  function  M(t)  are  given  by 

*(£)  =8  = r(*+l)(i||t||)-^(|t|) 
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andM(t  )=6  = T (i/+l  ) (±flt|)_l/ I„  (||t||)  . 

Proof : Any  real  vector  t can  be  written  as  Jt|  a where  a is 

a unit  vector  in  the  direction  of  t.  The  result 

g eii'u  = g en±h'u  = }ei»t|w  (i-w»)^dw 

J1  B(|>  ^+i) 

=r(i/+i)  (iitij-^j^iti) 

is  now  immediate  from  Result  3 and  formula  (8.4.10)  of 
Gradshteyn  and  Ryzhik(1980) . The  corresponding  result  for 
the  moment  generating  function  is  now  obvious  by  Result  2 
and  formula  (8.43.1). 

This  result  is  somewhat  more  significant  than  it 
appears.  If  we  partition  U into  Uj  and  U2  as  we  did  in 

Result  3,  it  should  be  clear  that  we  can  determine  the 
characteristic  function  and  moment  generating  function  of 
Hi  from  the  above  result.  If  we  let  be  a kxl  vector,  we 
can  write  a / as  a,=  (a/,0,)  where  0/  is  an  lxm-k  vector  of 
zeros,  thus  A„(|a|)  = g e-  ^ = g e ~ 1“1  It  should  be  clear 

that  ||a  II  = Ill'll- 


Result  6.  Let  Ui,U2,  . . . ,Un  be  independent  and 
distributed  random  vectors  with  Ux  — U ( Sm ) , and 


ident i cal ly 


let  Rn— 


t—i 


then 
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(i)  £ i R"U> 

t=l 

(ii)  U~U(Sm)  and 

(iii)  Rn  and  U are  independent, 

Proof : See  Mardia  (1979) , Chapter  15,  section  15.4. 

Result  7 . With  Rn  as  defined  in  Result  6,  we  have 
r2=r^+r 

n — m +2R  mRn-mCosS  where  Rm,  Rn  _ m and  cos#  are  mutually 
independent  and  t = cos#~Wi/ . 

Proof : Now 


7/t  Tl 

+ ^2  U,.  ~ RmU*  + Rn-mU**  + 2R„R„- mU'*U* * , 


i= 1 


t — m + 1 


where  Rm  is  the  length  of  y^U- , R„-m  the  length  of  ^ U< , U* 

* — 1 i=m+l 

and  U**  are  unit  vectors  . Hence 


R"=  (£n,)'(Ju<)  ~ Rn=Rm+R„_m  + 2 RmR„- 

i= 1 Xi=l 


uL*u, 


Since  U**  is  a unit  vector  then  the  density  of  t 
conditioning  on  U** , say  g (t  | U**)  , does  not  depend  on  U** 
and  by  Result  2 


(l-tV(m  3) 


-l<t<l 


«(*  III**) 


0 


e 1 sewhere . 


-14- 


It  follows  "that;  "the  unconditional  density  of  t is  the  same 
as  the  conditional  density,  that  is,  t~W^. 

We  conclude  this  chapter  by  summarizing  a paper  by 
Saw  (1986) . 


1 . 4 On  a Paper  by  Saw. 

Formula  (8.532)  in  Gradshteyn  and  Ryzhik(1980)  us  to 
write 

cr(coS  ») 

^m+n  k= 0 un 

where  R^J_|_n=Rm+R2+2RmRncos0 ; Rn_|_m  , Rn , and  Rm  are  the  sides 
of  a triangle  such  that  the  angle  between  Rm  and  R„  is  9, 
Cn(cos0)  is  Gegenbauer  polynomial. 

Again  in  Gradshteyn  and  Ryzhik(1980) , formula 
(7.323.1),  we  have  the  result  that 


Cn  (cos  9s)  sin 2U9  d#  = 


0 


0 


2~2t/7rf(2r+l) 


f(r+l): 


n=l ,2 , 


n=0. 


Noting  the  relationship  between  Aj/(t)  and  the  Bessel 


function  Jj/(t)  and  using  the  formulae  mentioned  above,  Saw 
arrived  at  the  apparently  simple  but  powerful  result 
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g Ai/(tR„)=  {A„(t)}“ 


(1.4.1) 


An  alternative  proof  of  (1.4.1)  is  to  let  U,  i=l,2,...,n 
be  independent  identically  distributed  random  unit  vectors 
then 


g (U1+U2+ . . . +U„)  _ g e t(W1+W2+.  . .+Wn) 


= {Ar(t)}". 


Now 


(u1+uJ+...+u„)  = |u,+u,+...tu.|.^;v;^> 


and  lUj+^f  . . . +Un||  is  independent  of  ^1,+T!+  ' ' '+^T?V  Using 

||iil+ii2+ • • • +LU 


the  result  of  Theorem  1 we  get 

« e"'(U‘+U2+'  ••+S«  = * = A1/(tRn)  . 

The  last  equality  is  immediate  on  noting  that  U~U(Sm) . 

The  usefulness  of  equation  (1.4.1)  lies  in  the  fact 
that  it  gives  a computational  method  for  determining  the 
distribution  of  Rn . We  now  develop  the  theory  as  given  by 
Saw . 

Let  Xf  be  a chi  random  variable  with  f degrees  of 
freed  om . It  is  straightforward  to  show  that 


.,2  k _ okr(k+|(f+2s)) 

Xf+2,  ~ * r(I(-f  + 2s))  ’ 


g 
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and  thus  it  should  be  apparent  that 


g 


A!,(txm+2a)  = 8 A, 

2 i 

= M (gin+s  , lm  : It2)  = 


(■tXnH-a,) 


V'  ('I-t2')  * r(^m)  r(k+l  (m+2s)  ) 

Q 2 k ! T(i (m+2s)  ) r (k+-|m) 


it2 

e2  M (-s  , |m  : -|t2)  =e9  M(-s,  P|:-0) 


where  0=lt2  and  M(a,b:z)  is  the  confluent  hypergeometric 
function.  The  fourth  equality  is  an  application  of  formula 
(9.212)  in  Gradshteyn  and  Eyzhik  (1980). 

A property  of  the  confluent  hypergeometric  function 
M(a,b:z)  is  that  when  a is  a negative  integer  or  when  a 
equals  zero  the  series  has  a finite  number  of  terms.  Thus 


g Ai/(txm)  — e^ . 1 
6 Ai/(tXm+2)=  e*.(l+fl/(lm)) 


This  can  be  written  in  matrix  form  as 


g Aj,(tXm) 
8 Al/(tXm  + 2) 


-e0A  D-1  6 = e6m, 
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where  A = 


10  0 
12  0. 
14  4 0 


D-1= 


1 

0 

0 


0 


i 

m 


0 

0 


0 1 \ 
m(m+2) 


9 


0_,=  ( 1 ,9,92,  . . . ) and  m = A D 19_ 


We  can  show  by  polynomial  division  that  A„(t.Jg)  can  be 

it2 

written  as  the  product  of  e2n  and  g(t2,m,n)  where  g(t2,m,n) 
is  a polynomial  in  t2 . Consequently  we  have 


g 


H2 


G(It2) 


5 


where  G(|t2)  = g(t2,m,n)|n.  We  can  express  the  polynomial 
G(|t2)  as  the  dot  product  of  two  vectors,  the  component 
of  the  first  vector  being  the  coefficient  associated  with 
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(^t2)*  and  "the  k1"'1  component  of  the  second  vector  being  the 
coefficient  associated  with  (It2)*,  thus  giving 


Using  the  fact  that  a/  = a/DA  *AD  1 we  obtain  the  result 
eda/9  = e'V  D A-1  A D_10  = eV  D A_1m, 


by  equation  (1.4.2). 

Thus  we  can  write  8 Aj/(tJ^  R„)  as  a linear  combination  of 

8 Ai/(tXro+2>)  s=0 ,1,2,  . . . 

Let  E be  the  infinite  diagonal  matrix 
E = diag  (1,2, 4, 8, . . .) 
and  p be  the  infinite  vector 

P=  (p(Xm<rc)  ,P(xm+2<rc)  , P(*m+4<rc)  , . . . ) , 
and  let  B be  the  matrix 


B = 


10  0 0 
-110  0 
1-210 


9 


then  A 1 


E 1 B,  giving  8 AI/(t^R„) 


e^a,DE-1Bm  . 


result 


From  the 
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g A„(tjf  R„)  = e^a'  D E-1  B m, 

we  conclude  that; 

P (Rn< r)  = a/D  E-1  B p. 

Thus  Saw  has  demonstrated  that  Tor  all  m>2  and  all  n>2 
we  can  write  the  distribution  Tunction  oT  R„  as  a linear 
combination  oT  the  distribution  Tunctions  oT  chi  random 
variables.  This  is  an  exact  relationship  involving  an 
inTinite  sum.  In  Chapter  II  we  will  mention  some  oT  the 

problems  one  encounters  in  evaluating  the  distribution 

Tunction  oT  Rn  using  the  method  described  above. 

The  useTulness  oT  equation  (1.4.1)  is  Turther  demonstrated 
by  considering  the  situation 

l 

Rn  = ||d1U1+d2U2+  . . .+dnUn|2 

where  d,-  i = l,2,...,n  are  constants  and  U,-  i = l,2,..,n  are 
independent  random  vectors  with  each  having  a uniTorm 
distribution  on  the  unit  m— sphere  . Then  Trom  equation 
(1.4.1)  we  have 

n 

g Aj,  (tcRn)  = IP  i/(ctd,)  (1.4.3) 

1 = 1 

and  by  "the  same  considerations  we  again  obtain  £ Al/(tcRn)  = 
9 t 

e a ,0  Tor  some  vector  a.  It  is  now  clear  that  we  can  write 
g Aj,(tcR„)  as  a linear  combination  oT  g Aj,  (txm_|_2jfc) 

k— 0 ,1,2,... 
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Saw  also  discussed  the  situation  in  which  the  d,- ’ s are 
independent  and  identically  distributed  with  emphasis  on  a 
Xm+2a  distribution.  Again  by  (1.4.1) 


S - {S  *m+2,} 

= { e%(-s,  r+1;  ^)}". 


From  our  prior  discussion  M(-s,j/+1,^)  is  a polynomial  of 
order  s in  9.  And  so,  our  polynomial  is  now  of  finite 

degree . 


CHAPTER  II 

A COMPUTATIONAL  STUDY 


2.1  On  a Paper  by  Saw  (Continued) 


The  result  that 

P (Rn<r)  = a'  D E_1  B p, 

gives  us  a method  Tor  computing  the  distribution  Tunction 
oT  Rn  n>2,  Tor  any  m>2 , by  using  a linear  combination  oT 
the  distribution  Tunctions  oT  some  chi  random  variables. 
Furthermore,  the  structure  oT  the  matrices  E-^  and  B,  and 
the  relationship  between  P(xy<r)  and  p(x/+2<r)  allow  us  to 
write  and  use  the  alternative  expression 


P(R„<r)  = a/D  d P(Xm<r) 


(2.1.1) 


where  d is  the  inTinite  vector,  d = (1  , d , <92 , . . . ) , and  d an 

operator  deTined  by  <9g(m)  = ||g(m+2)  - g(m)|. 

We  can  convey  an  idea  oT  the  ease  oT  the  use  oT 
equation  (2.1.1)  by  examining  how  to  obtain  the  vectors  and 
matrices  given  in  the  equality.  We  saw  that 


8 A„  (tcR„)  = jAj,  (tc)  = ^'n(^nm)  = e^a'fl. 


where  ^(0): 


00  u 

E ok- 

k=0 


k!  m(m+2) 


a_  1-4-2 

(m+2k-2)’ 


and 
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The  computation  of  the  coefficient  of  9k  in  is  best 

effected  by  repeated  use  of  an  algorithm  corresponding  to 
polynomial  multiplication.  We  can  consider  Vfr  ( ^ nm ) as  an 

infinite  vector  with  the  coefficient  of  0j_1  as  its  jth 
element  j=l,2,...  In  order  to  obtain  the  coefficient  of  9k 
in  the  expansion  of  we  select  the  first  k+1 

elements  of  the  vector  , and  on  making  n-1 

applications  of  the  algorithm,  we  obtain  the  first  k+1 
elements  of  the  vector  ^(^p3)  with  the  last  element  being 
the  coefficient  of  9k . The  power  series  for  e-^  is  known; 
we  can  now  obtain  the  elements  of  a/  by  polynomial 
multiplication  of  ^"(tt3)  by  e~ 9 and  letting  the 

coefficient  of  9 ^ k=0 , 1 , 2 , . . . in  the  resulting  polynomial 

be  the  k+1  element  of  a. 

The  diagonal  element,  dj.j.,  of  the  infinite  diagonal 

matrix  D is  known  and  can  be  written  as 

dfcfc=m(m+2) (m+4)  . . . (m+2k-2)  k=2,3,...  and  dn=l  . 

The  vectors  and  matrices  we  are  considering  each  have 
an  infinite  number  of  elements.  We  must  truncate  these 
vectors  and  matrices  to  q+1  terms,  q = 0,1,2,...;  Saw 

refers  to  the  value  obtained  as  the  q^^1  approximant  to 
P(Rn<r).  When  m=3 , the  density  function  of  Rn  is  known  in 
closed  form.  We  can  compare  the  values  of  the  various 
approximants  with  exact  distribution  of  R„  for  different 
values  of  r and  n.  This  is  done  in  Table  I,  which  is  taken 


from  Saw. 
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Result  2.1.1  is  exact . However,  in  order  to  use  "this 

result;  we  must;  truncate  "the  vectors  and  matrices  given  in 
the  expression.  Table  II  indicates  that  the  approximants 
do  not  rapidly  converge  to  the  exact  value  of  the 
distribution  function  of  R„  when  n is  small  (n<9) . 


2.2  A Computational  Study 


Examine  Table  II  of  section  2.1  and  note  that  for  m=10, 
n=2  and  r=0-5,  the  10^  approx  i mant  is  equal  to  -0.0001. 
We  obtain  negative  values  for  some  approximants  since  some 
of  the  coefficients,  in  the  linear  combination  of  the 
distribution  functions  of  chi  random  variables,  are 
negat i ve . 

A chi  random  variable  is  defined  on  the  positive  half 
of  the  real  line,  while  R„  is  defined  on  a finite  interval. 
The  slow  convergence  of  the  approximants  towards  the  exact 
value  indicates  that  we  must  include  a large  number  of 
terms  in  order  to  have  the  linear  combination  of  the 
distribution  functions  of  chi  random  variables  be  a good 
approximation  for  the  distribution  function  of  Rn . 

A part  of  our  research  was  to  obtain  a set  of  random 


oo 


var i abl es 


such  that 


k=l 


8 Aj,  (tR„)  — a0-{  Aj/  (t ) aj.6  Aj/(t  Y^.) 

k—  1 
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TABLE  I 


Tabulated  is 

diTf=|P{Rrj(  ^<r}  - P?{Rn^  ^<r}  | and  reldiff  = d iTT/P{R„  ^3^<r} 
Tor  some  specimen  values  oT  n,  r and  q. 


n = 3 

q = 

0 

diTT 

re 1 d i TT 

r = 0.5 

1 .0 

E-2 

4.8 

E-l 

r = 1 .0 

3.2 

E-2 

1.9 

E-l 

r = 1.5 

4.0 

E-2 

9.2 

E-2 

r = 2.0 

3.0 

E-2 

4.3 

E 

1 

II 

to 

Or 

1.7 

E-2 

1.8 

E-2 

r = 3). 

2.9 

E-2 

2.9 

E-2 

q = 10 

d i TT  re  Id iTT 

7.0  E-5  3.3  E-3 
6.9  E-4  4.2  E-3 
9.5  E-4  2.2  E-3 

7.4  E-4  1 .0  E-3 

1 .5  E-3  1.6  E-3 

2.1  E-3  2.1  E-3 


r 

r 

r 

r 

n 


r 

r 

r 

r 


4 q = 0 q = 5 

d iTT  re  Id iTT  d iTT  re  Id iTT 


0.5 

2.5 

E-3 

1 .4 

E-l 

4.9 

E-4 

2.7 

E-2 

1 .5 

3.5 

E-2 

1 . 1 

E-l 

9.8 

E-4 

3.0 

E-3 

2.5 

6.0 

E-3 

7.5 

E-3 

1 .9 

E-3 

2.4 

E-3 

3.5 

1.7 

E-2 

1 .8 

E-2 

2.3 

E-3 

2.3 

E-3 

6 

q = 
diTT 

0 

re 1 diTT 

a = 
diTT 

5 

re 1 d i TT 

0.5 

1.3 

E-3 

1.3 

E-l 

5.2 

E-5 

5.2 

E-3 

2.0 

2.2 

E-2 

5.4 

E-2 

4.3 

E-4 

1 . 1 

E-3 

3.5 

6.3 

E-3 

6.9 

E-3 

2.6 

E-4 

2.9 

E-4 

5.0 

4.5 

E-3 

4.5 

E-3 

1 .6 

E-4 

1 .6 

E-4 
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TABLE  II 


Comparison  ot  various  approximants  with  the  exact  value 


<N 

II 

C 

r = 

0.50 

dim 

= 

10 

20 

30 

40 

50 

Exact 

m = 

2 

0 . 1742 

0.1621 

0 . 1732 

0.1767 

1 . 1685 

0.1609 

m = 

3 

0.0681 

0.0644 

0.0667 

0.0633 

0.0621 

0.0625 

m = 

4 

0.0268 

0.0277 

0.0271 

0.0257 

0.0267 

0.0260 

m = 

5 

0.0104 

0.0120 

0.0112 

0.0117 

0.0110 

0.0112 

m = 

6 

0.0039 

0.0051 

0.0049 

0.0051 

0.0049 

0.0050 

m = 

7 

0.0013 

0.0021 

0.0022 

0.0022 

0 . 0023 

0.0022 

m = 

8 

0.0003 

0.0009 

0.0010 

0.0010 

0.0010 

0.0010 

m = 

9 

0.0000 

0.0004 

0.0005 

0.0004 

0.0005 

0.0005 

m = 

10 

-0.0001 

0.0002 

0 . 0002 

0.0002 

0.0002 

0.0002 

n = 2 

r = 

1 .00 

dim 

= 

10 

20 

30 

40 

50 

Exact 

m = 

2 

0.3224 

0 . 3384 

0 . 3386 

0.3294 

0 . 3337 

0 . 3333 

m = 

3 

0.2401 

0.2518 

0 . 2483 

0.2491 

0.2491 

0.2500 

m = 

4 

0 . 1915 

0.1989 

0 . 1940 

0 . 1942 

0 . 1964 

0 . 1955 

m = 

5 

0.1561 

0.1582 

0.1554 

0.1569 

0.1568 

0.1562 

m = 

6 

0 . 1283 

0 . 1267 

0.1257 

0 . 1270 

0 . 1268 

0.1266 

m = 

7 

0.1060 

0 . 1028 

0 . 1031 

0 . 1037 

0.1038 

0 . 1035 

m = 

8 

0.0880 

0 . 0844 

0.0853 

0.0854 

0 . 0851 

0 . 0852 

m = 

9 

0 . 0732 

0.0698 

0.0708 

0.0706 

0.0705 

0.0706 

m = 

10 

0.0612 

0.0581 

0.0589 

0.0586 

0.0586 

0.0587 
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Table  Il(continued) 
n=2  r=l . 4142 


d im= 
m = 2 
m = 3 
m = 4 
m = 5 
m = 6 
m = 7 
m = 8 
m = 9 
m = 10 


10 

0 . 4837 
0.4989 
0.5020 
0.5012 
0.5005 
0.5005 
0.5009 
0.5015 
0.5020 


20 

0.4769 
0.4928 
0.4922 
0.4951 
0.4972 
0.4977 
0.4979 
0.4983 
0 . 4987 


30 

0.4856 

0.4944 

0.4964 

0.4993 

0.4995 

0.4995 

0.4999 

0.5000 

0.4999 


40 

0.4934 

0.5006 

0.5012 

0.5013 

0.5007 

0.5005 

0.5005 

0.5005 

0.5002 


50 

0.4982 

0.5038 

0.5024 

0.5013 

0.5007 

0.5004 

0.5003 

0.5003 

0.5001 


Exact 

0.5000 

0.5000 

0.5000 

0.5000 

0.5000 

0.5000 

0.5000 

0.5000 

0.5000 
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oo 

where  aj.>  0 and  afc=l  • We  mush  report,  failure  in  this 

k= 1 

area  oF  the  research  but  we  now  outline  some  oF  our 
methods . 

For  n=2  and  any  rn>2  , it  is  known  that  has  a 

B^l(m-l) , l(m-l)j  distribution,  thus  an  intuitive  approach  is 
to  conduct  the  search  amongst  the  Family  oF  beta  random 
variables . 

Furthermore  ±(m-l)  = ±(m+j-3)  where  j=2  leads  us  to 

consider  the  set  oF  beta  random  variables 

oo 

(f (m-i) , |(m+j-3)|. 

j=2 

A second  approach  was  based  on  examining  the  First  two 
moments  oF  j2  . These  moments  can  be  obtained  by  a simple 
manipulation  oF  the  Function  Aj,(jf)  . Next  we  selected  a 
beta  random  variable;  its  First  two  moments  agree  with  the 
First  two  moments  oF  • Such  a random  variable  has  a 

b(| (mn-2) , | (mn-2) (n-l ) ) distribution.  Thus  a second  set  oF 
random  variables  For  which  we  could  investigate  the 
relationship  was 

{ B(  2ij(mj-2)  ’ j=2 ,3,4,  . . . 

The  method  oF  selecting  a beta  random  variable  by  matching 
the  First  two  moments  appears  in  a paper  by  N.L.  Johnson 
(1965) . 
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Table  III  gives  the  coefficient  of  (^b2)*,  k=0 , 1 , . . . , 6 , 

in  "the  expansion  of  when  v=\  and  n=2,3,...,8. 

Table  IV  gives  fhe  coefficient  of  (^b2)fc  k=0,l,...,6,  in 

the  expansion  of 

8 A„(t  B(  (jH_-2).(j-D)  ) 

where  v=\  and  j=2,3,4,  ...  ,8.  Note  fhaf  fhe  row 
corresponding  fo  n=2  and  fhaf  fhe  firsf  fwo  columns  of 
Table  III  are  fhe  firsf  fwo  columns  of  Table  IV. 

The  dafa  in  Table  III  and  Table  IV  can  be  considered  as 
fhe  elements  of  two  matrices,  each  having  seven  rows  and 
seven  columns.  We  can  now  obtain  the  coefficients  ELk 
k=0 ,1,2,. ..,6,  for  a f ixed  val ue  of  n by  post mu It i ply i ng 
the  inverse  of  the  transpose  of  the  matrix  in  Table  IV  by 
the  transpose  the  appropriate  row  from  Table  III.  The 
coefficients  obtained  from  solving  this  linear  system  were 
negative  at  times;  and  so  we  did  not  continue  with  this 
method  of  examining  the  moments  of  sets  of  beta  random 
variables  . 


2 . 3 Cone 1 us i on 


Our 
above . 


investigations  were  more  extensive  than  indicated 
The  program  which  gave  the  output  for  the  data  of 


Table  III  and  Table  IV  with  a few  minor  changes 
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can  be  used  "to  examine  higher  dimensions  and  different:  sets 
of  Beta  random  variables.  This  was  done.  In  higher 
dimensions  fhe  results  were  more  appealing  but  coefficients 
(a^  , k=0 , 1 , 2 , . . ) , were  at  times,  negative  and  their 
absolute  values  could  be  larger  than  unity. 

Our  failure  to  find  a suitable  set  of  random  variables 
does  not  indicate  that  the  problem  does  not  have  a solution 
but  rather  that  we  must  develop  a better  methodology.  Our 
methods  were  not  much  more  than  educated  guesswork;  perhaps 
there  are  functional  relationships  we  could  have  exploited 
but  which  we  were  not  aware  of.  The  fact  that  the  linear 
combination  of  distribution  functions  of  chi  random 
variables  can  be  replaced  by  a single  and  easily  computed 
distribution  function,  when  n=2 , leads  us  to  believe  that 
further  investigations  of  this  problem  may  be  profitable. 
A next  method  of  attack  is  needed  in  order  to  make  further 
progress  towards  the  solution  of  this  problem. 
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TABLE  III 


Coeft ic  ient 


in  the  expansion  ot  -[Aj, 


m=3 


k=2 


k=3 


k=4 


n=2 

n=3 

n=4 

n=5 

n=6 

n-7 

n=8 


1 . 0000000000  0 . 3333333333 
1.0000000000  0.2222222222 
1.0000000000  0.1666666666 
1.0000000000  0.1333333333 
1.0000000000  0.1111111111 
1.0000000000  0.0952380952 
1 . 0000000000  0 . 8333333333 


0 . 0444444444 
0.0213991769 
0.0125000000 
0 . 0081777777 
0.0057613168 
0.0042759961 
0.0032986111 


k=6 


k=7 


k=8 


n=2 

n=3 

n=4 

n=5 

n=6 

n=7 

n=8 


0.0001410934 

0.0000432893 

0.0000170855 

0.0000079770 

0.0000041934 

0.0000024059 

0.0000014758 


0.0000042755 

0.0000011099 

0.0000003757 

0.0000001525 

0.0000000707 

0.0000000361 

0.0000000200 


0.00000000939 

0.00000000211 

0.00000000062 

0.00000000022 

0.00000000009 

0.00000000004 

0.00000000002 


k=5 

1.0031746031 

>.0011902947 

>.0005621693 

>.0003074708 

>.0001206350 

>.0001206350 

>.0000826719 
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TABLE  IV 
2 

Coefficient  (%-)*  in  the  expansion  of 
g { A„  (t 

2J  2 J ' J 

m=3 


k=2 


k=3 


k=4 


k=5 


CN 

II 

1.0000000000 

0 . 3333333333 

j=3 

1 .0000000000 

0.2222222222 

II 

•n 

1 .0000000000 

0.1666666666 

j =5 

1 .0000000000 

0.1333333333 

j=6 

1.0000000000 

0.1111111111 

j=7 

1 .0000000000 

0.0952380880 

j=8 

1 .0000000000 

0 . 0833333333 

0.0444444444  0.0031746031 
0.0213991771  0.0011734036 
0.0125000000  0.0005527210 
0.0081777778  0.0003023716 
0.0057613173  0.0001828989 
0.0042759955  0.0001188830 
0.0032986111  0.0000815590 


k=6 


k=7 


k=8 


II 

to 

0.0001410934 

0.0000042755 

j =3 

0.0000417878 

0.0000010468 

C_i. 

II 

0.0000163129 

0.0000003460 

j =5 

0.0000076034 

0 . 0000001395 

co 

li 

•n 

0 . 0000040028 

0.0000000646 

j=7 

0.0000023021 

0.0000000332 

0.00000000939 

0.00000001947 

0.00000000055 

0.00000000019 

0.00000000008 

0.00000000003 


j =8  0.0000014159  0.0000000184  0.00000000001 


CHAPTER  III 

ON  THE  DISTRIBUTION  OF  R„ 
3 . 1 Literature  Review 


Determination  of  the  probability  density  of  Rn , the 
length  of  the  sum  of  n independent  random  vectors,  is  a 
problem  which  has  been  investigated  by  a considerable 
number  of  workers.  Rayleigh  (1905)  gave  an  approximate 
solution  in  the  distribution  of  the  resultant  when  each  of 
the  independent  random  vectors  has  a uniform  distribution 
on  the  circle.  At  this  time  his  work  was  motivated  by  a 
question  posed  by  Pearson  (1905).  At  a later  period, 
Rayleigh  gave  an  exact  solution  to  the  three  dimensional 
case.  Then  his  investigations  were  more  extensive  as  he 
also  examined  the  probability  density  of  R„  in  one  and  two 
d imensions . 

Kluvyer  (1905)  gave  an  exact  solution  to  the  two 
dimensional  case.  However,  his  answer  is  in  the  form  of  an 
integral,  with  the  integral  being  a product  of  Bessel 
functions.  Evaluation  of  this  integral,  even  with  the  use 
of  numerical  techniques,  is  difficult.  G.N.  Watson  (1944), 
on  page  421,  in  his  treatise  on  the  Theory  of  Bessel 
Functions  extends  Kluvyer ’s  method  to  an  arbitrary  number 
of  dimensions;  Watson’s  solution  is  given  by 
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P r | eL^  9 • . . 9 cln)  — 


, , ''T  im.,  n f J1  m_,(a<t)  1 

{r(»}"  X{  (ir-fc)2  Jim(rt)n|  — jdt  (3.1.1) 

0 ( n-a.-'t ) 


(la.t) 


where  ai : i = l,2,...,n  is  the  length  of  the  i*3*1  random 

vector . 

Like  Kluvyer’s  solution  the  right  side  of  equation 
(3.1.1)  is  seen  to  be  an  integral  which  involves  the 
product  of  Bessel  functions.  Even  when  computed  by  the  use 
of  numerical  integrations  methods,  convergence  of  this 
integral  is  slow;  thus,  an  alternative  method  of  computing 
the  distribution  function  of  Rn  is  of  some  importance  to 
workers  in  this  area. 

In  this  chapter  we  develop  such  a method.  We  first 
introduce  the  situations  where  an  alternative  method  is 
already  known  in  section  3-2  and  section  3-3.  In  section 
3.2  we  give  the  distribution  of  R2  for  any  m>2  whereas  in 
section  (3-3)  we  discuss  the  distribution  of  Rn  (m=3) . In 
section  (3.4)  we  give  the  distribution  of  R„  for  odd  values 
of  m . 

3 . 2 On  the  Distribution  of  R„ 


In  Chapter  I we  mentioned  the  distribution  of  arU,  where 
a is  a unit  vector  and  U has  a uniform  distribution  on  the 


m-sphere . 
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It;  is  this  result  which  allows  us  ho  derive  "the  probability 
density  function  of  R2  for  m>2  . 

Theorem 1.  Let  Hi  > and  U2  be  independent  identically 

distributed  random  vectors  with  U^UCSm).  Let  ax , and  a2  be 
positive  real  numbers  and  R2  be  the  resultant  of  a^fa^ ; 
then  the  probability  density  function  of  R2  is  given  by 


h(r)  = 


aia2B(i  , \ (m-1)  ) 


0 


[ 


(r^~(al  ~a-2)2)((a-l  + a-2)2~  r2)-i  2*-m  3) 


2 0^2 


] 

| ai - a2  | < r <(a!+a2) 
e 1 sewhere . 


Proof : The  scalar  product  (aJJ!  + a2U2),(a1U1  + a2U2)  gives 


Rl 


2 2 

aj^  + a2  + 2a1a2cos0, 


where  9 is  the  angle  between  Ux  and  U2 . We  saw  that 
T = cos  9~VI1' . 

The  quantity  a2  + a2  + 2a1a2cos0  attains  its 
value  when  0=0  and  its  minimum  value  when  0=7r . 

Thus  for  | a1-a2  | < r <(a!+a2) 


maximum 


P(R2  <r)  = P (R2  < r2)  = Pfcos0  < ), 


2 2 2 
r -aj-aj 

2 ax  a2 


and  on  denoting 


by  m(r)  we  get 
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m(r) 


P^cos^  < m(r)  J 


(1-x2/ 


km-3) 


dx . 


Denoting  "the  distribution  function  of  R2  by  H(r)  we  write 
0 r<  | a1-a2  | 

H(r)=|  P^cosfl  <m(r)  ) |a!-a2|<r<  (a1+a2) 

1 r > (ai+a2)  . 


As  a result  of  the  Liebenitz  theorem  for 
differentiating  under  the  integral  sign,  we  obtain  the 
density  of  R2  as 


a-ia2  B(i,i(m-1))  lai-a2l<  r < (ai+a2) 

h(r)  = J 

0 el sewhere , 


and  on  noting  that  the  expression  in  square  brackets  is  a 
difference  of  two  squares  we  write  h(r)  as 


h(r)  = 


aia2B(I,i(m-1)) 


r(r2—  (ai~a2)  2)  ( (a1+a2)2-r2)  |(m~3) 

1 2a:a2  ‘ 


Letting  ax  = a2  = 1 gives  density  function  of  the  length 
of  the  resultant  of  the  sum  of  two  independent  random 
vectors  on  the  unit  sphere  as 


h(r)  = B(l,l(-1  ))Sr’<4-r‘» 


(m-3) 


0 < r < 2 . 


We  now  develop  the  three  dimensional  case. 
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3 . 3 The  Three  Dimensional  Problem 


In  section  1-4  we  saw  "that; 


8 A1/('tRn)  — |Aj,  (t ) (3.3.1) 

where  v = ^m-1  and  m denotes  the  dimensionality  of  the 
problem.  We  can  obtain  a similar  result  with  the  function 
Vj,(t),  where  Vj/(t)  = A„(it),thus 

S V„(tRn)  = {Vj,(t)}n  (3.3.2) 


Let  X be  a random  variable  having  a uniform 
distribution  on  the  interval  [-1,1]  then  the  characteristic 
function  of  X,  g eltx,  is  given  by 


g 


i tx 


eit-e~it 
2 it 


Vi(t) 


(3.3.3) 


Equation  (3.3.2)  and  equation  (3.3.3)  are  useful  in 
deriving  the  probability  density  and  probability 
distribution  function  of  Rn  for  the  three  dimensional  case. 


Theorem  2 . Let  Y be  a random  variable  with  probability 
density  function,  g(y) , satisfying 


(i)  g(y)  = g(-y) > 
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>0 

o 

V 

>> 

(ii)  g(1)  (y ) = | 

=0 

o 

II 

>> 

\ 

<0 

y < 0. 

Let  <j>y(  t)  = S 

iby 

e be 

the  characteristic 

function 

of  Y. 

Then  -2rg^^(r) 

i s the 

probability  density 

function 

of  a 

positive  random 

variable 

say  R,  such  that 

S Vl(tR)  = <t>y(  t)  . 
2 


Proof : Clearly  -2rg^'\r)  > 0,  and 


oo  oo 

j -2rg(1)(r)dr  = -2r  g(r)  + 2jg(r)dr  = 1, 

0 0 


hence  -2rg(1)(r)  , r>0 , is  a probability 

Also 


g Vi(tr) 
2 


OO  . 

'eltr-e-ltr 
2 i t r 

0 


!-2rgW(r)]dr 


density  function. 


oo 


/„it  „-itr\  fi-v 

) *‘V>  dr 


0 


/0itr  ^-itr\ 

= ( S 5t ) SC")  I 


OO 


oo 

0 


(eitr-  e-itr)  g(r)  dr 


0 


OO 


eitr  g(r) 


-OO 


dr  = <t>Y  (t)  , 


which  is  the  result  we  wanted  to  establish. 
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Coro  1 1 ary . If  we  denote  "the  probability  distribution 

function  of  R by  H(r)  and  that  of  Y by  G(y)  then 

H(r)  = 2G(r)  - 1 - 2rg(r) . 

Proof : We  know  that 


H(r)  = P(R  < r) 


g(1)(t)  dt 


0 


= -2t  g(t)  \rQ  + 2 

0 

= -2r  g(r)  + 2 G(r) 

= 2G ( r)  - 1 - 2rg ( r) 


g(t)  dt 
- 2G (0) 


which  is  the  required  result. 


By  letting  G(r)  = 1 - G(r)  and  H(r)  = 1 - H(r)  we  get 

the  alternative  expression 


|H(r)  = G(r)  + r g(r)  . 


Let  Y = Y1 X,  where  Xi  i = l , 2 , . . . ,n  are  independent 
i=l 

identically  distributed  random  variables  with  the 

probability  distribution  function  of  Xx  given  by 
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5 -l<Xl<l 

f(xi)  = | 

^ el sewhe  re . 

Then,  from  the  discussion  preceding  Theorem  2,  we  have  the 
result 


tfyCO  = {vi(t)}n  = 8 ViCtRn)  . 

k 2 J 2 

Theorem  2 and  the  corollary  following  this  theorem  give  us 
the  distribution  of  Rn  for  the  three  dimensional  case  we 
have 

h (r)  = -2r  g(1)(r)  r>0  , 

and 

iH(r)  = G(r)  + r g(r)  , 

where  h(r)  and  H(r)  are  the  probability  density  function 
and  probability  distribution  function  of  R„ . The  functions 
g(r)  and  G(r)  are  the  density  function  and  distribution 
respectively,  of  the  random  variable  Y. 

We  conclude  this  section  by  giving  the  density  and 
distribution  function  of  Y.  A derivation  can  be  found  in 
the  works  of  Feller  (1971),  Cramer  (1946)  and  Rao  (1973) 
amongst  others. 

The  density  function  g(y)  is  given  by 


1 

2"  ( n- 1 ) ! 


t (-l)J'(n)[y+n-2jr"1 

i= o J 


g(y)  = 


(3.3.4) 
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and  distribution  function  by 


G(y)  = E (-1)J(n)[y+n-2j]n  (3.3.5) 

^ n ' i=o  j ' 


where  -n  < y < n and 


0 (y+n-2 j ) < 0 

[y+n-2  j]=j 

(y+n-2 j ) (y+n  — 2 j ) > 0. 

3.4  Higher  Dimension  (m,  odd). 

The  literature  on  the  subject  of  Bessel  functions  is 
quite  rich.  These  functions  have  been  studied  in  depth  and 
their  properties  well  tabulated;  see  for  example  Watson 
(1944).  We  mentioned  the  relationship  between  V„(t)  and 
the  Bessel  function  J„(t)  in  Chapter  I.  We  now  show  that 
Vj/(f)  like  Jj/(t)  , satisfies  a differential  difference 
equation.  It  is  this  property  of  Vjy(t)  which  allows  us  to 
develop  the  results  of  this  section. 


Lemma.  The  function  V^(t)  satisfies  the  differential 

difference  equation  v2uVu  (tw)  \=2vvi21'~1V  u_  1 (tw)  (3.4.1) 


Proof . The  series  representation  of  w2I/Vl/(tw)  is 


w-v„(tw)  =|0 


2 k ,2r+2ifc 

w 
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Since  "the  series  is  uniformly  convergent;  we  can  interchange 
the  operations  of  summation  and  differentiation  to  obtain 


d / 2ry  si  ,\  1 r (r+1 ) (2r+2k)  . . o / j,  % 

dw  tW  V^Ctw)|  - Eo  Ic ! r (r+k+l ) (-1) 


2k 


W 


2J'+2ifc-l 


which  on  simplifying  gives 


,2V+2k-l 


= 2i/ 


w2"_1V 


(tw)  , 


which  is  the  result  we  wanted  to  prove. 


Multiplying  both  sides  of  the  equality  in  equation 
(3 .4.1)  by  w~^  we  get 

( = 2rw2I/-2VI/_1(tw)  (3.4.2) 

which  leads  us  to  conclude  that 

2*r(r-l)  . . . (r-k+l)w2I/-2*  V^.^tw)  (3.4.3) 

Let  W be  a random  variable  with  a symmetric  probability 
density  functions,  say  g(w),  then  the  characteristic 
function  of  W can  be  written  as 

OO  00  . . 

S eltW  = J e 1 twg  (w)  dw  = J (- ^ )g(w)  dw 

-oo  0 
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00 

f V j (tw)  g(w)  dw. 

1 2 


Furthermore 
identical ly 


if  W=^W,-,  where  the  W, ( 
«— 1 

distributed  random  variables 


s are  independent 
with  W^W1',  then 


S e 


i tW 


oo 


0 


V j (tw) g (w) dw 
"2 


= [V„(t)]n, 


and  on  denoting  the  density  of  R„  by  h(r),  we  have  the 
result 

op  oo 

|Vt/  (tr)h(r)dr  = |vi/(t)|"  = f V i (tw)  g (w)  dw  (3.4.4) 

0 o' 

Equations  (3.4.1)  through  (3.4.4)  will  be  useful  in 
determining  h(r),  the  probability  density  function  of  R„ , 

n 

from  a knowledge  of  the  density  of  wibh  W.^W^. 

i = l 

In  this  section  the  number  of  dimension  m is  related  to 
the  variable  k by  the  equality  m = 2k+l  = 3,5,7, . . . 

We  now  develop  our  main  result  by  the  principle  of 
mathematical  induction. 


Theorem  3.  Let  gk  (x)  k=l  ,2,  . . . be  a sequence  of  functions 
such  that  the  first  k derivatives  of  gj.  (x)  exist  and  let 


the  functions 
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V!  (W-0  (tX) [fek)ig*  (X> 11 


j=0,l ,2, . . . ,k 


be  integrable , "then 


oo 


V 1(tx)gl(x)dx= 


-OO 


(-1)' 


1.3... (2k-l) 


00 

1 VK-I)(tx)x2‘[(^)*g‘(x)ldx 


(3.4.5) 


Proof:  Observe  that 


oo  oo 

[ V|(tx)  gi(x)  dx  = | ^[x  Vi(tx)]  gx  (x)  dx 


-OO 


oo 

= X Vi(tx)  gx(x)  I | X Vi(tx)  g(11)(x)dx 


-oo 


oo 


Vi(-tx)  x2[(^g^)gi(x)]  dx. 


Thus  "the  result  is  "true  for  k=l  . (Note  that  the  function 
xg[1}  (x)  is  similar  to  the  density  of  R„  in  section  3.3.) 
Assume  that  the  result  is  true  for  k=i , that  is 


(-I)'' 


oo 


V,(tx)g,(x)dx  =1.3,5v;.V(2i_1)  j I(^;),gi(x)]dx 


-OO 


It;  is  clear  from  equation  (3.4.1)  that 


35<fi+1  vj(!,+1)(tx))  = (2i+1>  x’‘  vj(*-0(tj0 
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thus  ■for  k=  i + 1 we  can  write 


oo 


v i(tx)g,+1(x)dx= 
'2 


(-1)’ 


-OO 


1.3...  (2 i - 1 )J 


-oo 


(-i),+i  f 

1.3.5...  (2 i + 1 )J  Vi(2i+1)^tx)  x 

-oo 


2i  + 2 


1.3.5. . . (2 i+1 ) 


— KSU 


which  is  the  result  we  set  out  to  prove;  thus,  by  the 
principle  of  mathematical  induction,  the  result  is  true  for 
all  positive  integers  k. 

Coro 1 1 ary ; Let  gfc(x)  be  a sequence  of  even  functions  that 
satisfies  the  conditions  of  Theorem  3,  then 

oo  oo 


is  an  even  function 


since  it  is  the  product  of  even  functions. 
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The  second  integral  of  equation  (3.4.5)  may  now  be  written 
as 


(-i)' 


oo 


1 .3.\  . (2k-l)  Vi(2t-i)(tx)  dx 

-OO 

k 00 

1.3.  ^ . (2k- 1 ) IVi(2Jfc-i)^tX^  x2i[(^)fcg^x^  dx 


which  completes  the  proof. 


Coro  1 1 ary : Let  Wj^W^  i = 1 ,2,  . . . ,n  be  independent  random 

variables.  Let  g(w)  be  the  probability  density  function  of 

n 

and  h(w)  denote  the  probability  density  function  of  R„ 

i—i 

in  the  m-d imens ional  case,  then 


h(w) 


2(-l) 

1.3.5. 


\ (m-0 
. . (m-2) 


Proof : First  define  the  even  function  h*(w)  by 


h*(w)=h*(-w)= 


hO) 


w>o  . 


From  the  discussion  preceding  Theorem  2,  we  know  that 


op  oo  oo 

I V 1 (tw)h*(w)dw  — f Vj  (tw)h(w)dw  = f V t 

J 2m  — 1 J 2m_1  J “2 

-oo  q — oo 


V 1(tw)g(w)dw, 
~2 


As  a result  of  Theorem  3,  we  write  this  equality  as 
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30  OO 

I V-I<*»)h*(w)dw»J  V, 


r/_d_N5(m-0 

L V.  ,Ax. i) 


-OO 


-OO 


(m-2) L v wdw 


g(w)] dw 


which  can  be  written  as 


J V-I(tw)  [h*(w)-1.3.5...(m-2)^)2(m  g(w)]dw=0  (3.4.6) 


We  know  "that:  the  integrand  in  equation  (3.4.6)  is  a 

continuous  function  in  the  variables  t and  w.  Multiplying 
this  equation  by  t2l/  we  can  then  differentiate  both  sides 
of  this  equality  with  respect  to  t.  Equation  (3.4.3)  is 
indicative  of  how  we  must  proceed  from  here.  We  get 


° = 


OO 


-OO 


*2VVi  (tw)(h*(w)  - 

jm-i  V 

(-W2) 


2\2^m  0 
(m- 

2N5(m-0 


1 .3.  . . (m-2)  LVwdw)2<^  ^S(w)])  dw- 


/ 2\2v  / 1/  \ 

Since  the  functions  h*(w)  and  ^ — ? _/ (_d  ) 2 /w\  1 

v 7 1 .3.  . . (m-2)l>wdw'  '/ 

are  even  functions,  it  is  possible  to  write  the  equality 
above  as 


oo 


I 

0 


eitw  [h*(w) 


/ 2 \ 2 ( m 0 1 /•  \ 


which  leads  us  to  conclude  that 


h*(w)  = 


(_w2)K-  0 

1.3... (m-2) 


for  all  values  of  w. 
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In  particular  tor  w > 0 we  get 


/ 1/  , \ 

= 2 few.?.(-3){(5a5),("'l)«oo} 


(3.4.7) 


which  is  the  result  we  set  out  to  prove 


In  our  proof  of  the  corollaries  following  Theorem  3,  we 
used  the  result  that  g(w)  is  an  even  function  in  order  to 
show  that 


h(w)  — 2 


(-w2) 


2x5(m-1) 


1.3.  . . (m-2) l vwdw 


-2){^ 


7 ' 


g 


(w)} 


rn  >3 


This  is  important.  The  function  V x (tw)  is  an  ev 

~2 

function;  therefore,  when  g(w)  is  an  even  function  we  get 


(tw)g(w)dw  = 


oo 


0 


1 (tw)g(w) 
"2 


dw 


oo 


= 2 


cos(tw)  g(w)dw, 


0 


that  is,  we  can  consider  this  integral  as  a Fourier  cosine 
transformation.  Let  g(t)  be  the  Fourier  cosine  transform 
o~f  2g(w)  then  Tranter  (1956)  gives  the  inversion  formula  as 


oo 

2g(w)  = Jg(t)  cos  wt  dt  . 

0 


If  g(w)  is  an  odd  an  integrable  function,  then 
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oo 

K 


V_1(-tw)  g(w)  dw  = 0; 
2 


it  follows  "that  "for  an  arbitrary  function  g(w)  "the  function 
obtained  from  the  inversion  need  not  be  unique. 

S i nee 


1 — im 


vim = Gtw)  2 Ji  ,N  r (|m) 

O"*  * rj  m 1 


we  have  that 


oo 


0 


V‘(tw)  dw 


oo 

■J 

o 


1 2\  \ ( m~  0 l /•  x 


OO  - 

= (it)1^mr(»[j1  (tw)  w f(-iy 

J om  ^ -L  • O • • • 

0 


Km_i)  i 


1-1  rr  d N (im-0 


(5T 2)W  g(w)]dw. 


. i 
1~2m 


Ignoring  (|t)  2 , the  above  integral  by  equation  (1.5)  of 

Tranter  (1956)  is  the  Hankel  transform  of 


1.3. . . (m-2) 


If  we  denote  the  Hankel  transform  of  this  function  by  f (t) , 
then  the  inversion  formula  is  given  as 
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oo 


f(t)  t Jx  (tw ) dt 


0 


l-Jm 


Multiplying  f(t)  by  r(im)(it)  2 and  equating  this 
expression  to  8 Vx  (tRn)  , we  get 

sm  — 1 


r (§m)  (it)  _5mf(t)  = {Vi^^t)}^: 


r (im) J i , (t) 


m — 1 


(^)2 


i 

r)  YVb  1 


-r 


so  that 


f(t)  = 


_ V2 


uJm_1  r(§m)  Jl  ,(t) 


(it) 


rO) 


C- 


(it)2 


m — 1 


thus 


i °9  i . Jj  (t) 

h(w)  = w^CrCjm))"-1!  (Jt)sm"lt  J,  (tw)[J^ ]»  dt 

J 2m_1  ,,  . =m-l 

0 (it)2 


(it): 


This  is  equation  (2.2.21)  of  Watson  (1982).  Thus  we  have 
the  result  that 


2(-w2) 


2n  \ (m_1) 


1.3.  . . (m-2) L vwdw 


re  d nK™-1) 


g(w)] 


1 ^ ^ j ^ (^  ) 

=w2m(r(im))"-1[(it)2m_1t  Jj  (tw)[J^ ]"  dt 

J «~‘  (it)^1 


0 


We  therefore  have  an  alternative  expression  for  the 
probability  distribution  of  R„  for  odd  values  of  m (m>3)  . 

We  now  give  the  probability  distribution  function  of  R„ 
for  the  2m+l  dimensional  case. 
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Coro 1 1 ary : With  g(w)  an  even  function  satisfying  the 

conditions  of  Theorem  3,  then 


1 


3 


(-I)"  f i 2m  / d 

5.  . . (2m-l)J  ^tdt 

-w 


)mg(t)dt 


m — 1 


= 2G(w)  - 1 + £(-l) 
:= o 


m+j^j-m+l 


r(l> 


( d n 

f(m- j+|)  vwdw; 


m — j — l 


g(w)  t 


w 

where  G(w)  = | g(w)  dw . 

- w 

Proof : Consider  “the  integral 

w 

|t2„(_d_rg(t)  dt>  (3  4 8) 

-w 

on  integrating  by  parts  we  get 

w 

(3.4.9) 

-w 

Ignoring  the  constant  (2m-l),  then  the  integral  in  equation 
(3.4.9)  is  that  of  equation  (3.4.8)  with  m replaced  by  (m- 
1) . Integrating  by  parts  again  we  get 

2w2m'1(^)m_1S(w)  - 2(2m-l)w2"-3(^)"-’g(w)  + 

w 

2(2m-l)(2m-3)J  t2— 4 (^)  — 2g(t)  dt . 

-w 
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Repeat;  the  operation  of  integrating  by  parts  until  the 
integral 


( - 1 ) m_1 (2m- 1 ) (2m-3)  . . .3 


w 

tg(1)(t)  dt 


-w 


is  obtained.  Integrating  by  parts  we  get 


(- 1 ) m_1 (2m- 1 ) (2tn-3)  . . .3  [2wg(w) 


w 


g(w)  dw] 


= (-1 ) m_1  (2m- 1 ) (2m-3)  . . .3  [2wg(w)  +1-2  G(w)]. 


We  can  express  (2m-l) (2m-3) . . . (2m-2j+l)  as 

2Jf  (m  +|) 

f(m  - j + I)' 


This  allows  us  t.o  write 


w 

' t2m(tft)m  g^)dt  = 

-w 


m — 1 

E 

j= o 


f(m-j+i)  vwdw 


m — j — 1 


g(w)  + (-l)m-12m 


r(m+|) 

r(i) 


[1 -2G (w) ] . 


It  follows  that 


(-l)m 

1.3.5... (2m- 1 ) 


W 

-w 
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m — 1 

2 G(w)  - 1 + £ (-1) 

• — n 


3 = 0 


which  is  the  result  we  wanted  to  prove. 

Denoting  the  distribution  function  of  R„  by  H(r)  and 
using  the  result  that 


For  completeness  we  give  a second  method  of  determining 
h (w)  from  a knowledge  of  g(w)  . This  method  is  based  on  a 
result  by  Watson.  We  do  not  reproduce  Watson’s  proof  here; 
we  merely  state  the  result  and  show  that  the  two  functions 
satisfy  equation  (3.4.T) . In  Chapter  IV  we  state  and  prove 
a general  result,  which  specifically  includes  Watson’s 
result . 


r 


we  get 


H (r)  = 2G(r)  - 1 


Letting  H = l-H(r)  and  G(r)  = l-H(r),  we  get 


In  equation  (2.15),  Watson  (1982)  gives  the  result 
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n 


un(y)  = |g(r)h(^)4f, 


y 


where  un(y)  is  "the  sum  of  n independent  W17  random 
variables,  g(r)  is  the  density  of  the  resultant  R„ , and 


unit  vector  and  U~U(Sm). 

Expressing  Watson’s  result  in  our  notations  we  write 


where  the  number  of  dimension  m is  given  by  m=2k+l . 
Equation  (3.4.6)  indicates  how  we  are  to  show  the 
relationship  between  h (w)  and  g(w)  . Differentiating  both 
sides  of  equation  (3.4.10)  with  respect  to  y and  then 
dividing  by  -y  we  get 


On  repeating  the  operations  of  first  differentiating 
and  then  dividing  by  -y  a further  (k-2)  time,  we  get 


h (t ) is  the  probability  density  of  aAj,  where  a is  a fixed 


n 


(3.4.10) 


Now,  by  using  a theorem  of  calculus,  on  differentiating 
with  respect  to  y we  get 
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-^-V-Wvi  = F(^k+^2  h(y) 

^ ' dyl  ydy/  sCy;  p(i)  v2*-i> 


which  can  be  written  as 


h(y) 


2 (- 1 ) * 

1.3.5... (2k- 1 ) 


)*  s(y) ; 


"this  is  equat  ion  (3. 4. 7), in  which  w is  re  pi  aced  by  y . 

In  chapter  IV  we  express  equation  (3.4.7)  in  a form 
that  is  suitable  tor  computational  purposes. 


3 . 5 Some  Problems 


It  is  clear  that  tor  higher  odd  dimensions  (m>3)  , we 
now  have  an  alternative  method  tor  computing  the 
probability  density  function  and  the  probability 
distribution  function  of  the  length  R„  of  the  sum  of 
independent  vectors  on  the  unit  m — sphere . How  does  this 
method  compare  with  Watson’s  general  result? 

The  Fourier  inversion  theorem  and  the  convolution 
theorem  are  powerful  methods  in  determining  the  probability 
density  function  of  the  sum  of  independent  random 
variables.  However,  when  the  random  variables  are  defined 


on  a finite  range,  the  probability  density  function  of  the 
sum  consists  of  several  components,  each  of  which  is  valid 
for  a specific  subinterval  of  the  range. 
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The  mathematics  required  to  evaluate  the  probability 

n 

density  function  of  ^^W,-  by  the  convolution  method  is 

i= 1 

rather  straightforward.  This  exercise  requires  the 

integration  of  a polynomial,  the  identification  of  the 
subintervals,  and  the  determination  of  the  limits  of 
i ntegrat i on . 

Since  the  probability  density  function  of  a higher 
value  of  n is  computed  from  that  of  the  next  lower  value, 
great  care  must  be  taken  to  avoid  an  error  in  calculation, 
since  a mistake  at  an  earlier  stage  can  invalidate  much  of 
the  subsequent  work.  The  degree  of  the  polynomial  which  we 
must  integrate  increases  as  n increases  and  also  as  the 
number  of  dimensions  increase.  The  number  of  subintervals 
increases  as  the  values  of  n increase. 

Watson’s  result  of  section  3-1  is  quite  general;  here 
the  vectors  are  of  equal  length.  In  Chapter  IV  we  will 
present  our  result;  which  will  also  encompass  the  situation 
where  the  random  vectors  are  of  fixed  but  unequal  length. 


CHAPTER  IV 

SPHERICAL  DISTRIBUTIONS 
4 . 1 Introduct ion 


Most  of  the  literature  on  multivariate  statistical 
analysis  is  devoted  to  the  theory  of  the  multivariate 
normal  distribution.  One  reason  that  the  multivariate 
Gaussian  distribution  plays  such  a central  role  in 
statistical  analysis  is  that  many  of  the  results  that  are 
valid  for  the  univariate  normal  distribution  can  be  easily 
extended  to  the  multivariate  case  as  for  example,  the 
univariate  central  limit  theorem.  As  in  the  case  of  the 

univariate  normal  distribution,  if  the  statistician 

believes  that  the  underlying  random  process  is  not 
adequately  described  by  his  multivariate  Gaussian  model, 
he  must  then  formulate  alternative  distributions. 

One  approach  is  to  formulate  a general  class  of 
distributions  that  contain  the  multivariate  normal 
distribution  as  a special  case.  The  advantage  of  this 
approach  is  that  many  of  the  results  which  have  been 
established  for  the  normal  distribution  will  hold  true  for 
members  of  this  general  class;  thus,  many  of  the 

statistical  tests  will  be  valid,  even  though  the  assumption 
of  a normal  distribution  has  not  been  satisfied. 


-56- 


-57- 


One  class  of  density  functions  which  has  received  some 
attention  is  the  family  of  spherical  densities.  Examples 
of  spherical  densities  are 

(1)  the  multivariate  normal  distribution,  X~N(0,I), 

(2)  the  uniform  distribution  on  the  unit  sphere  U~U(Sm) , 

(3)  the  uniform  distribution  on  the  unit  m-ball, 

(4)  the  multivariate  t distribution  with  n degrees  of 
freedom  and  density  function 


ffe)  = -r(bd-)  (1  + „ 

f (in)  (nir)  '2,m 


and 


(5)  the  multivariate  Cauchy  distribution  (the  multivariate 
t-distribut ion  with  n=l) . 

In  section  4.2  we  give  some  results  for  a random  vector 
X having  a spherical  density.  In  sections  4.3  and  4.4,  we 
give  the  relationship  between  the  probability  density 
function  of  the  length  of  an  m-component  spherical  vector 
and  the  probability  density  function  of  a single  element  of 
that  vector.  In  sections  4.5,  4.6,  and  section  4.7,  we 
give  some  examples  using  the  multivariate  normal 
distribution,  the  multivariate  Cauchy  distribution,  the 
uniform  distribution  on  the  m-sphere , and  the  uniform 


distribution  of  the  m-ball. 
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4.2  Some  Distributional  Properties 


In  "this  section  we  give  some  of  the  properties  of  a 
random  vector  having  a spherical  density.  It  is  these 
properties  which  we  use  to  develop  the  results  in  section 
4 . 3 and  4.4. 

In  section  1.3,  we  constructed  a probability  space  by 
lett i ng 

Cm  (C  >t)={u : u e Sm , C'xj.t} 
and  letting 

A =(t  |cm(C,t):  C e Sm,  t c [-1,1]} 

be  the  smallest  (7-algebra  generated  by  the  collection  of 
all  ’’caps”.  Define  the  non-negative  measure  U such  that 


0 


1 <t  <oo 


c|cm(C,  t ) } = 1 -A (t ) = 


(l-x2) 

B(^,r+l) 


dx 


-1 <t <1 


1 elsewhere. 

A random  vector  X is  said  to  have  the  property  of 
spherical  symmetry  if  P(C'x  < t)  is  independent  of  C when  C 
£ Sm . An  immediate  result  is 


Theorem — 1.  Let  X be  a random  vector  which  is  spherically 


symmetric,  then 


-59- 


(i)  pq~U(Sm)  and 


(ii) 

Proof 


■p^pand  ||X||  are 

c^x 

Let  P(^fxf^w’ 


stochastically  independent. 
lX|<r}=H (w , r)  and  !{.}  be  the 


i nd i cator 


function  of  the  event  {.}.  Since 
H(w,r)=Sugx  I{|^<w}.I{IXS<r} 
=«X  A(w).I{|X|<r} 

=A(w)P(|X|<r) 


the  result  follows. 


Theorem  1 gives  a representation  of  a spherically 
symmetric  random  vector. 

Corollary  . If  we  let  X be  a spherically  symmetric  random 
vector  then  we  have  the  following  distributional  result 
X ~R  U, 

where  R - |X|  > U~U(Sm)  , and  R and  U are  stochastically 

i ndependent . 
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Some  writers  describe  a,  spherically  symmetric  random 
vector  by  one  of  the  following  definitions: 

Definition  1 . If  the  random  vectors  X and  PX  have  the 
same  density  function  for  all  fixed  orthogonal  matrices  P, 
then  the  vector  X is  said  to  have  a spherical  distribution. 

Definition  2.  If  the  probability  density  function  of  a 
random  vector  X can  be  written  as 

f (x)  = g(x'x)  , 

then  the  vector  X has  a spherical  distribution. 

That  these  two  definitions  are  equivalent  should  be 
obvious  when  considering  Theorem  1 and  Corollary  2,  however 
we  give  a proof  of  the  result. 

Theorem  2.  Definition  1 and  Definition  2 are  equivalent. 

Proof : Let  the  probability  density  function  of  X be  given 

by 

f (x)  = g (x'x)  . 

Make  the  transformation  Y=PX , where  P is  a fixed  orthogonal 
matrix  such  that  | P | =1 . Since  the  Jacobian  of  this 

transformation  is  unity,  we  obtain  the  density  function 
h (y)  of  Y as 
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h(y)  = g(y'P'Py)  = g(y'y) ; 

"thus,  X and  Y=PX  have  -the  same  distribution  for  every  fixed 
orthogonal  matrix  P. 

Suppose  X and  PX  have  the  same  probability  density 
function  for  every  fixed  orthogonal  matrix  P.  Denote  the 
density  function  of  X by  m(x)  . Make  the  transformation 
Y=PX  with  P a fixed  orthogonal  matrix.  Since  X and  PX  have 
the  same  density  function,  it  follows  that  m(x)  satisfies 
the  functional  equation 

m(x)  = m(P'x)  (4.2.1) 

for  every  fixed  orthogonal  matrix  P. 

IT  2^—0.  then  equation  (4.2.1)  is  true;  hence,  we  need 
only  examine  the  case  of  x^O . 

Every  real  vector  x can  be  written  as  ||x||a,  where  a is  a 
unit  vector  in  the  direction  of  x-  Thus  the  functional 
equation  (4.2.1)  is  written  as 

m(  |x|a)  = m ( ||x|  P'a)  . 

With  x fixed,  P;  can  be  a fixed  orthogonal  matrix,  with  a 
as  its  first  row,  giving 


m(  | x || a)  = m(  pxlej) 
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where  is  a vector  with  unity  as  its  first  element  and 

all  other  elements  equal  to  zero.  Thus,  for  every  fixed 
vector  we  have  x 

m(  | x 1 a)  = m(  |x|£1)=g4,  (|x|)  , 

giving  m (x)  = g*  (x'x)  for  every  fixed  real  vector  x . We 

can  now  conclude  that  the  two  definitions  are  equivalents. 

These  two  characterization  of  a random  vector  x having 
a spherical  distribution  are  useful  in  proving  some  of  the 
distributional  properties  of  this  family  of  distributions. 

Theorem  3.  Let  X.i>X2,  . . . ,X„  be  independent,  identically 
distributed  random  vectors,  with  Xi  having  a spherical 
distribution,  then  X.1+X.2+  • • • +Xn  has  a spherical 
distribution . 

Proof : For  every  fixed  orthogonal  matrix  P,  we  have  that 

P (X1+X2+  • • •+2£n)=  PX.1+PX.2+  . . . +PX„ 

= X1+X2+  • • • +Yn 
~Xl+X2+  . . . +Xn  , 
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where  Yf  i=l,2,...,n  are  independent,  identically 

distributed  random  vectors  with  Y^ ~ X,  • This  completes  the 
proof . 

Theorem 4.  Partition  the  spherical  random  vector  X1  into 

j X/2)  , where  X^  is  a ixmt  vector  and  X^  is  a lxm-mj 
vector,  then  X\  has  a spherical  distribution. 

Proof : Our  proof  is  complete  if  we  can  show  that  for  every 

fixed  orthogonal  rr^x  mj  matrix  Q,  then  Q)£i  ~Xj . Consider  the 
matrix  P,  where 

Q Oi 

p= 

0/  I 

where  Q is  any  orthogonal  matrix;  0X  is  an  m:  x m-mt  matrix 
whose  elements  are  all  zero  and  I is  the  (m-m1)x(m-m1) 
identity  matrix.  P is  an  orthogonal  matrix;  thus,  PX~X 
gives  Q Xi~Xx  for  every  fixed  orthogonal  matrix  Q. 

These  two  properties  are  useful,  the  first  tells  us 
that  the  sum  of  n independent  random  vectors,  having  a 
uniform  distribution  on  the  m— sphere,  has  a spherical 
distribution.  Any  result  we  can  develop  for  this  class  of 
distributions  is  applicable  to  the  sum  of  independent 
random  vectors  with  a uniform  distribution  on  the  unit 


sphere . 
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The  second  result,  Theorem  4,  tells  us  that  if  X has  a 
spherical  distribution  then  any  subvector  has  a spherical 
distribution  as  well.  This  property  gives  us  the  results 
in  section  4.3. 

It  is  interesting  to  note  that  we  could  have  obtained  the 
proof  of  Theorem  3 and  Theorem  4 by  other  methods  as  for 
example,  by  the  use  of  characteristic  functions.  We  now 
give  some  further  results  on  spherical  distributions. 

Theorem  5.  Let  X be  the  random  vector  with  probability 
density  function 


f(x)  = cmh(x'x) 


and  let  U = and  R = |X|,  then 


(i)  R and  U are  independent, 

(ii)  U~U(Sm)  and 

(iii)  the  density  of  R is  given  by 


g(r) 


r(|m) 


Jh(r) 


r>0  . 


Proof : In 

X~N  (0  , I)  . 


Chapter  I we  proved  this  result  for  the  case  of 
The  proof  here  is  identical  except  that  instead 


of  e 


we  now  have  the  function  h (x?x)  . 
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Kelker  (1970)  partitioned  X'  into  (X'x , x'2)  , with  Xx  a 

9X1  component  vector . He  then  uses  the  characteristic 
function  and  a polar  transformation  to  obtain  the  joint 
density  of  and  r2=|X2|  as 


f Oi,r2) 


cm 


l 

2m  q — 1 m — q — 1 , 2 

ri  r2  g(r^+rj) 

r(sci)r(im-ici) 


The  next  result  is  taken  from  Lord  (1954) . 


Theorem  6 . 
distribution 


Let  X be  a random  vector  having  a spherical 
then  the  characteristic  function  is  given  by 


<K/0  = g e - 


1 

Tyfn 

r Jim_1(rvOs(r)dr> 

2 

where  p=|t| , r=|x|| , and  g(|xfl)  is  the  probability  density 

function  of  X. 


l 

2m 


= (2*) 2 p 


-2m+i 


Proof ; Make  the  transformation  U 


Mx|. 


The 


Jacobian  of  this  transformation 
characteristic  function  is  given  by 


Rm_1 


g e — = grgu 


«-ie  i rt  u 


l s 


The 
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Transform  U to  U*,  where  PU  = U*and  P is  an 

matrix  with  its  first  column  being  i , to  get 

t 


grgu  rm-1e  * r— *— 


Srg  * rm-1e 1 r^Ul 
r u,* 


where  p=||t|| . The  characteristic  function  of  U*  is 
given  by 


> *e 
Ul 


i rpuf  , , l-irn 

= 2 (l^)  2 Ji  (rp)T^m) 
2 


Using  result:  (iii)  oT  Theorem  5 we  get: 


i—  77m 


KP)=  Sr  2(|r p)  2 Jim_1(r/,)r(Im) 


1 1 1 

. 2m  l-~m  f 

(27r)  P r J j (r/?)g(r)dr , 

J 2m_1 

0 


which  is  the  result  we  wanted  to  prove 


Lord  points  out  that  <^(/>)  is  a function  of  p 
same  process  will  give  the  inversion  formula 


f(x)  = S ( r)  = (2tt  ) 


-jm+i 


OO 


0 


1 

ry  TTi 

P Ji  ,(rP)<t>(p)dp , 
2m~ 1 


orthogonal 


known  and 


only;  the 
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which  as  Lord  states  is  the  probability  d 
m variables.  However  we  can  obtain  the 
length  of  X,  by  using  a transformation . 


ensity  f 
density 


unction  of 
of  R,  the 


Theorem  7.  Let  X be  a random  vector  with  a spherical 
distribution,  then  the  probability  density  function  of  R=||X|| 
is  given  by 

1 00  i 

h ( r)  = 2^m+1{r(lm)}-1  j (rp) ( rp)  0 (p)  dp 

o 2m 

Proof : From  Theorem  6 we  have  the  result 


1 1 OO  1 

^(ti»t2 t m)  = <j>(p)  = (2ir)  5 p 5 [ r5  Jj  (rp)g(r)dr, 

J 2m-1 

0 

Then  by  the  Fourier  inversion  theorem  we  obtain  f(x)  as, 


f(-)  = (i)mj  el~“  *(£)  dtx . . . dt, 

Rm 


Make  the  transformation  U and  R = |X|| , then 


f(r’  = (i)mrm  1 j eir--^(i)  dtx  . . . dtm . 

Rm 

Transform  t to  ( p,a ),  where  «=Aj  and  p=|t|,  then 


00  • / 

f(r’  — ) = (^r)m  r”*  ^ je  P ~-pm~l  <j){p)  dad p, 

0 Sm 
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Lst  P be  an  orthogonal  matrix  whose  first  column  is  u_  then 
transform  a to  g*=Pg ; thus,  the  integral 

f ipra/u 

J e --da 

Sm 

becomes 


je1"™*  da*. 

Sm 

We  can  regard  this  integral  as  the  characteristic  function 
of  a single  element  of  a random  vector  that  has  a uniform 
distribution  on  the  unit  m— sphere , giving 


ipra^ 


1 

2 m 


da*  = 2tt  2 (Ipr) 


1—  4m 

Jlm  l(r^ 
2m  1 


>m 


The  marginal  density  of  r is  obtained  by  integrating  the 
joint  density  f(r,  u)with  respect  to  the  variables  of  u. 
We  have 


fO,  u) 


l 00  i 


m — 1 


0 


Ji  , ( rP ) •<£(»  dp. 

jm-l 


Upon  integrating  with  respect  to  the  variable  u,  we  get 


OO 


f(r)  = 


nfn-1 


r(|m)  o 


~m 

P Ji  , (rp)<j)(p)dp . 

2m  — 1 


This  is  the  inversion  formula  for  the  probability  density 
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function  of  fhe  random  variable  R,  "the  length  of  a 
spherically  symmetric  random  vector. 


In  Chapter  III  we  mentioned  Watson’s  result  without 
giving  a proof  for  it.  The  following  is  a more  general 
statement  of  Watson’s  result  and  develop  Watson’s  result  as 
a coro 1 1 ary  . 


Theorem  8 . 


Let  T be  a positive  random  variable  such  that 


T~{b ( jk , i(m-k)) |2 , and  let  T be  independent  of  R,  then  the 
distribution  of  Z = RT  is  given  by 


oo 


s(z)  =- 


ifc-i 
2 2 


rak)  0 


\k 


where  *(/>)  = S e1^*  and  R is  the  length  of  X,  a spherically 
symmetric  random  vector. 


Proof : Make  the  transformation  Z=RT  and  S=T-2.  The 

Jacobian  of  this  transformation  is  seen  to  be  |S_1 . The 

regions  0<R<oo  and  0<T<1  are  transformed  into  the  regions 
0<Z<oo  and  l<S<oo.  To  obtain  the  marginal  density  of  Z we 
must  integrate  with  respect  to  s. 

Thus  , 


oo  oo 


g(z)=- 


I! 


1 1 
-zm  Ttm 

P z2 


(2  r (Im)  B (Ik  ? Ifm-lcj)!  o 


'im+*(s-l)*(m"‘"° 


• JIm_!  ds. 


-70- 


Interchanging  the  order  of  integration  and  letting  r=lm-l 
we  get 


g(z)  = 


(2t/r(r+l)B(ik,^+l-ik))J) 


OO  oo  x 


, „ V-\k 

(s-1)  Jj/(/9zfs)ds  dp 


But  the  result  that 


°?  — —1/ 

jx  2 (x-l/  1 Jj/  (a-\|x)  dx  = r(/i)2Aia-I/JJ/_/i(a) 

1 

can  be  found  in  several  texts,  such  as  Watson’s  (1944) 
Letting  p — v— |k+l  and  a - pz  we  get 


S(z)~  l ._} [ P2  z2  <j)(p)  J,  (pz)dp. 

22  f(ik)  0 2 


the  result  we  wanted  to  establish. 


In  Chapter  I , we  saw  that  o/U  has  the  same  distribution 
for  every  unit  vector  a when  U has  a uniform  distribution 
on  the  unit  sphere.  With  this  fact  in  mind  it  should  be 
clear  that  the  characteristic  function  of  any  subvector  of 
U has  the  same  form  as  the  characteristic  function  of  JJ. 
We  know  that  for  a spherical  random  vector  X,  then  X ~ RU 
where  R=|X|  and  U~U(Sm).  It  is  obvious  that  the 
characteristic  function  of  any  subvector  of  X has  the  same 
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form  as  the  characteristic  function  of  X-  For  example,  if 
X~N(0,I),  then 

„ i-t'x  -4  t't 

8 e1--  = e 2 

for  all  t e Sm . 

Theorem  8 tells  us  that  if  we  multiply  the  length  of  an 
m-component  vector  having  a spherical  distribution  by  a 
{b (gk , g (m-k) random  variable,  then  the  resulting  random 
quantity  has  the  same  distribution  as  the  length  of  any  r- 
component  subvector  of  the  original  vector. 


Coro 1 1 ary  . Denote  the  length  of  the  m-component  vector  by 
R and  its  probability  density  function  by  h(r).  Denote  the 
length  of  any  k-components  of  this  vector  by  Rt  and  its 
density  by  gj.(r)  then  gj.(r)  and  h(r)  satisfy  the  integral 
equat i on 


8*0) 


oo 


h(v) 

— 2 


/ 2 2 \ 2 ( m k — 2') 

(v  -r  ) 

B(ik,i(m-k)) 


dv 


(4.2.2) 


Proof : In  Theorem  8 we  saw  that  Rjfe~RT  where 

T ~ {B(ik,i(m-k)}2  and  l<k<m-l  and  is  independent  of  R. 

Make  the  transformation  Y=RT  and  V=R  and  then  obtain  the 
marginal  distribution  of  Y by  integrating  the 
probability  density  function  with  respect  to  v. 


j o int 
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Letting  k=l  in  equation  (4.2.2)  we  get 


SiO) 


oo 


h(v)  ( 

m — 2 


2 2 \ 2 ( m *0 
v -r2)** 


dv, 


which  is  the  equation  found  in  Watson  (1982) . 

As  mentioned  in  our  introduction,  many  of  the  results 
which  can  be  proved  for  a random  vector  having  a 
multivariate  normal  distribution  can  also  be  proved  for  a 
random  vector  having  a spherical  distribution.  We  now 
state  and  prove  two  such  results.  These  results  are  taken 
from  Mu i rhead (1982) . 


Theorem 9.  Let  X have  an  m— variate  spherical  distribution 

with  P (X  = 0) =0 . 


(i)  If  W— where  a e Sm  and  a1  a=l  , then 


V = ^(m-l) 


W 


J(i-w2) 


has  a t — d i st r i but i on  with  m— 1 degrees  of  freedom. 

(ii)  If  B is  an  mxm  symmetric  idempotent  matrix  of  rank  k 
then  S=  |x||2  Bas  a beta  distribution  with  parameters  |k  and 
i(m-k) . 


Proof : Observe  that  in  results  (i)  and  (ii)  the  random 

variables  V and  S are  functions  of  the  random  vector 

IX! 

which  we  know  has  a uniform  distribution  on  the  unit  m- 
sphere . 
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Without  loss  of  generality  we  can  consider* 
a multivariate  normal  distribution,  Z~N(0,I) 
(1 ,0, . . . ,0) . Th  en  we  can  write 


iz|  uhere  2 

and  take  a ' 


has 

as 


V =■ 


- (m-1)  Zi 


l 

2\  2 


(m 

E z?) 

i=2 


V is  the  ratio  of  a univariate  N(0,1)  random  variable  and 
the  square  root  of  an  independent  Xm-i/(m-l)  random 
variable;  therefore,  V has  a tm_j  distribution. 

To  prove  result  (ii),  note  that 


- xi 
|xf  |X| 


and  so  we  can  again  assume  that 


S 


_Z_ 

»Z| 


where  Z has  a multivariate  normal  distribution  Z~N(0,I). 
Since  B is  a symmetric  idempotent  matrix,  we  can  find  an 
orthogonal  matrix  H such  that 


HBH= 


I* 

O' 


0 

0 


giving 
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S ' «i|B  III  - III  H'  H B H'  H B H'  H Ill- 


Making  the  transformation  Y=HZ  we  get 


K _ 

E Y? 
s = 

EY,2 


i=l 


It  is  well  known  that  if  Xa  and  are  independent  then 


xl 


Xa  + AJ 


which  completes  the  proof. 

This  next  result  was  given  by  Kariya  and  Eaton  (1977)  . 
Here  1=  ( 1 , 1 , . . . , 1 ) and  {l}={kl:  k e Ft1}  ; thus,  {jl}  is  the 

span  of  1.  The  proof  was  taken  from  Muirhead  (1982). 

Theorem — 10.  Let  X and  Y be  two  independent  random  vectors 

where  X has  an  m-variate  spherical  distribution  with 
— £0  — 0 and  Y has  any  distribution  whatsoever,  with 
P(Y£l)=0.  If 


r = 


E (X4-X)  (Y,-Y) 

i=i 


E (X.'-X)2  E(Y,-Y)2 
*=1  »=1 


m X • m Y • 1 

where  X = £ and  Y = E HT  then  (m-2)2  T r 

<=1  i=1  -J(l-r2) 

t m-2  distribution,  here  m>2. 


has  the 


-75- 


Proof  ; Let  M — i i l1;  "then  r can  be  written  as 


r = 


X'(I-M)Y 


J[X'(I-M)X  Y'(I-M)Y]‘ 


Since  I-M  is  idempotent  of  rank  n-1  , there  is  an  orthogonal 
matrix  H such  that 


H ( I -M) H= 


I„-i 


O' 


0 


0 


where  I„-i  is  an  identity 


matrix . 

Put  U=HX  and  V=HY  and  partition  U and  V as 


U = 


. j:  = [ 


where  U*  and  V*  are  m-lxl  column  vectors,  then 


r = 


U'H(I-M)H,V 


J[U,H(I-M)H,U  V'H(I-M)h'V] 


U*'V* 


12  * V* 


1^*1  l¥*l 


Since  X has  a spherical  distribution,  we  know  that  U=HX  has 
a spherical  distribution  and  any  subvector  U*  of  U al 
has  a spherical  distribution. 


so 
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Fu^thermore  H*  is  independent  of  V* . Conditioning  on  V* 
and  by  ( i)  of  Theorem  11.  with  a=i^.  then 


(m-2) 


1 

2 


r 

J(l-r2) 


has  the  tm_2  distribution  and  the  proof  is  complete. 


Since  for  any  spherically  symmetric  random  vector  X, 
X~RU  such  that  R and  U are  independently  distributed  with 
U~U(Sm) , it  is  clear  that  any  function  satisfying 

k (^O  =k ( i'll)  for  0<r <00  satisfies  k(X)~k(U)  . Examine  Theorem 

9 and  Theorem  10  and  observe  that  replacing  X,  by  RX 

leaves  the  variables  W and  S of  Theorem  9 and  R of  Theorem 

10  unchanged. 


4 . 3 A Distributional  Property 

Let  Y be  a random  variable  with  probability  density 
function  g(y ) satisfying 

( i)  g(y)  = g(-y)>0  and 
(i  i)  g(y)  = 0 y>l  , 

then 


g 


i ty 
e ^ 


oc 

8 E 


*=0 


(it)* 

k\ 


k= 0 


(it)2* 

2k! 


f*2k 


where  /J2Jfe=g  Y2* . 
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Knowing  that;  Y~W^  , we  can  easily  show  "that; 


^2  k 


(2k- 1 ) (2k-3) . . .3.1 
(m+2k-2)  . . . m 


k = 1,2,... 


and  hence  obtain  the  characteristic  function  of  the  random 
variable  Y as 


S e ( it)  . 


With  X a random  vector  having  a spherical  distribution 
and  with  a a unit  vector,  we  have 


g e^— = g 7— 


S tfm(itR) 


oo 


tfm(itR)dHm(R) 

0 

where  Hm(R)  is  the  probability  distribution  function  of  the 
random  variable  R=|X|.  But  the  distribution  of  n'u  does  not 
depend  on  the  unit  vector  a;  therefore,  a may  be  a unit 
vector  with  each  of  its  first  k elements  not  equal  to  zero 
and  the  last  (m-k)  elements  are  all  equal  to  zero. 
Partition  X into  (X^jX^)  where  X*i  is  a lxk  row  vector  and 
X'2  is  a lxm-k  row  vector  and  partition  a'  into  (a'i  ,a'2)  , 
then 


oo 

Jtfm(itR)dHm(R)  = g 

0 
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g eitRjba,1U1 


where  Rj.— flXjj,  and  is  a lxk  unit;  vector  having  a uniform 

distribution  on  the  k-sphere  and  is  independent  of  Rk . 


independent  of  Rt . 


This  gives 


oo 


0 

where  Hi.(r)  is  the  probability  distribution  function  of  the 
length  of  any  k components  of  the  vector  X. 

For  any  m-component  random  vector  having  a spherical 
distribution  we  have  the  result  that 


for  all  values  of  t. 


Consider  the  function 


this  function  with  respect  to  r we  observe  that 


thus,  on  integrating  by  parts  we  get 
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oo 

^2m+l)J  ^2m+i(ir"fc)h2m+i(r)dr=(2rn+l)  f r2m*2m+1  ( i rt ) • r“ 2mh2m+1  ( r ) d 
0 0 


OO 

= -J  r2m+1^2m+3  ( i rt)  ^{r-2mh2m+1  (r)  }dr 
0 


for  all  values  of  t.  It  follows  that 


h2m+3(r)  (2m+l) 


r2m+i  d rr-2mh 
dr  L 


2 m-M 


(r» 


“ Gm+l)  h2m+l  ( r ) 


r h2m+l(r) 

(2m+l ) 


(4.3.2) 


This  equality  allows  us  to  express  the  density  function  of 
the  length  of  a spherical  random  vector  in  terms  of  the 
density  functions  of  the  absolute  value  of  a single  element 
of  that  vector , if  the  number  of  dimensions  is  odd. 


Theorem  1 1 . Let  X be  a (2m+l)  dimensional  random  vector 
having  a spherical  distribution.  Denote  the  density 
function  of  the  length  of  any  k>l  components  of  the  vector 
X by  ht(r)  then 


m — 1 


O)  = i £ 


j= 0 


/m-L  (2m  2-j)!  >+1 . 0+0, 

C j ' (2m-l)!  (~2r>  h (r) 


(4.3.3) 


where  h0)(r)  = (^)Jh(r)  , and  h(r)  is  the  probability 

density  function  of  a single  component. 
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Proof  : Letting  m=l  in  equation  (4.3.3)  we  have 


h3(r)  = - r h(1)(r)  . 


Assume  “the  result  is  true  tor  m=s . For  m=s+l  we  have 


h’<+3<»  = (sfl)  h3.+l<»  - (srn)  r C+.O)- 


We  write  r h"’+1(r)  as  rh“'+,(r)  = ! [ ^ !h(l)(r) 


+ 


V1  (S  1)  • (2s-2-k)  ! . 2,i  0 -v  ifc+l,(*  + l)r  w 

jj=i  k ! (s-k)  ! (2s- 1)  !^-2^  (k+k-2s)r  h (r)  + 


(S— 1)  (S  1)  • ( n\(  »+l  l (*+l)  ( a -1 

U-lJ(2s-l)  !(_2)  r h 


and  2sh2j+1(r)  as 


‘ H V)g^  • -»rh(1,(r)  tf2s(s-)  £» (_ar)  - 


.h(t+1)(r)] 


On  substituting  for  2sh2j+1(r)  and  for  rh^+1(r),  we  write 


2sh!<+1(r)  - r h<',>+1(r)  =J  [ - 2 rh  “>  ( r) 


(s-1) 


kT(s-k)  ! (2s-l)  ; (2s-2-k)  ! (-2)- * (k2+ k-2s)  ]r 


i+1h(t+1)(r)  - 


(s-1) 

(2s-l) 


i(-2) 


»r«+i 


h(>+1)(r)] 


(4.3.4a) 
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The  expression  in  the  inner  square  brackets  reduces  to 


(2s-k)  \(S\  , 0xJ;+l 


2s 


m-2>' 


On  dividing  this  quantity  through  by  (2s+l)  we  get 


(2s-k)  !/s\  *+i 

(2s+l)  ilk)^)  • 


Using  this  result  to  simplify  equation  (4.3.4a)  we  get 

WO=i£(“)  h“+1>(r)  (4.3.4) 

The  result  is  true  for  m=s+l  and  hence  it  is  true  for  all 
positive  integers  m by  the  principle  of  mathematical 
induction . 


Lemma : Let  a*.  _ ^ ^ )^~(2m+2)  ^ ‘ bhen  we  have 


(k+l)a‘+‘  + a*  = (k!t)'  (2II0  ! ' (-2>‘ 


Proof : The  proof  is  straightforward  and  has  been  omitted, 


Theorem  12.  Let  H2m+3  (y)  =Jh2m+3  (t ) dt , then 

0 


H - mvtVrn+1\(2m+2-k)  ' ,(*- 1),  . 

2m+3(y)  E(  k J ('2m+2)  ! ( 2)  y h (y)  (4.3.5) 
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Proof  : 
form 


Examine  equation  (4.3.4),  then  H2m+3 (y)  is  of  the 


H2m+3(y)  = \ 5>*y*h(t  °(y) 

A;— 0 


where 


h°(y)  = h (y)  and  h ^ *\y)  =J  h(t)dt. 

0 

Differentiating  H2m+3(y)  we  get 

H2m+3  (y)  = h2m+3(y)=i[mE1okafcy*-1h(t~l)(y)  + a*y‘h  (0  (y)] 

= ^ £o((k+1)a,+i^)yth(°y  + am+iym+lh^m+1\y)]  • 

In  Theorem  11,  we  established  that 


h2m+3  (y  ) 


(2m-k) ! 
(2m+l) ! 


(-2y) 


*+1h 


0+0,  x 

(y) ; 


thus,  equating  the  coefficient  of  ym+1 , we  get 


am+l 


m ! 

(2m+l) ! 


(-2) 


m+l 


(m+1) ! 
(2m+2) ! 


2(-2)m+1 


Using  the  lemma  preceding  the  theorem  we  determine  the 
value  of  auk  as 


a* 


C?1) 

k 


(2m+2-k)_! 
(2m+2) ! 


(-2) 


k 


k=m , m-1 , ,1 
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On  subst  ifut  ing , we  see  -that  ax=-2  and  thus  a0=2  since  ax  + 
a0  = 0 . Thus , 


H2m+3  (y  ) = H(y) 


\(2m+2-k) ! 
/ (2m+2) ! 


( -2y)  fch  (*_l)  (y)  , 


which  is  equat  ion  (4.3.5). 


4.4.  Solving  the  Integral  Equation 


In  section  4.2  we  demonstrated  that  the  probability 
density  function  of  the  length  of  an  m-component  spherical 
vector  is  related  to  the  probability  density  function  of 
the  absolute  value  of  a single  element  of  that  vector  by  an 
integral  equation.  We  are  able  to  solve  this  equation 
both  for  even  and  odd  values  of  m.  The  solution  we 
obtained  for  when  m is  odd  is  identical  to  that  given  in 
section  4.3.  We  first  solve  this  equation  for  odd  values 
of  m . 


Theorem  13 . 
equat i on 

S2ife+i(r) 


Let  h(r)  and 

DO 

f h(y)  (y2-r2)t~1 

J.  y2i-1B  (1 , k) 


Safc+i(r> 

dy 


satisfy  the  integral 


(4.4.1) 


where  m=2k+l , then 


1 *^Vk-l\(2k-2-j)! 

2 jk<S  j J (2k- 1 ) ! 


(-2r)'+1  g0+l)(y) 


h (r) 
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Proo~f : Letting  k=l  , the  result 

d i ff  e rent  i at  i ng  both  sides  of  equation 
to  r.  Assume  that  the  result  is  true 
k=s+l  we  have 


is  immediate  on 
(4.4.1)  with  respect 
for  k=s . Now  for 


&2S+3 


OO 

« ■ B 


(y)(y2-r2)- 
y2,+1B(l,s+l) 


dy 


Differentiating  both  sides  of  this  equation  with  respect  to 
r we  get 


oo 


D fv'l  = fh(y)  (y2-r2) *-1 
5+3  k ; J y2s+i  B(l,s+1)  ' 


-2sr  dy, 


with  a little  arithmetic  this  equality  can  be  written 


as 


oo 

(rifl)  C2sS2.+3(r)-rgD)+3(r>]  = J }J  ' dy 

f y D \2 5 s/ 

Denoting  (3^)  [2sga,+8(r) -rg^+8(r)  ] by  g*(r)  we  get 


(r)  = J 


h(y ) (y2- r2)  * 1 
y2*-1B(| , s) 


dy , 


but  by  the  induction  process 


h(r)  = 


1 ,'r'1  (2s- j -2)  ! 

2 jk 0 (2s-l)  ! 


("j1)  (-2r)>+‘g„<'+'>(r) 


(4.4.2) 


The  Liebenitz  rule  for  repeated  differentiation  of  a 
product  allows  us  to  write 
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S*('+l)(r)=(-2^7j[(2s-j-l)g<',+‘>(r)  - r g'^O)]  ; 


on  subst itut i ng  for  g*^J  ^ (r)  in  equation  (4.4.2),  we  get 


h ("  r')  — 1 V (~^s  j • (s-li,-  q„\1+i  Vfn  • (i+i)  f \ 

^ ) 2(2s+l)  (2s-l)  ! V j 2 ^ [(2s-j-l)g2s+3  (r) 


-«*g2+?(r)]  . 


is  easy  to  show  that  the  coefficient  of 
j =0 ,l,...,s-l  in  the  above  express ion  is 

(rin)(i)(p  gSfi<-2)i+' 

while  that  of  r,+1  is 


(2s+l) 


l(s)  s!(-2) 


8-\~l 


giving  the  final  equality  as 


h (r)  = l 


= * yg,  (liilj!  <2s-J>!  (j)  (-2y)i+14,.«,(r). 


By  the  principle  of  mathematical  induction,  the  result  is 
true  for  all  positive  integers  k. 


We  now  investigate 
values  of  m.  Although 
difficult  than  when  m is 


this  integral  equation  for 
this  situation  is  slightly 
odd,  any  elementary  text  on 


even 


more 
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i ntegral  equations  will  illustrate  how  to  solve  this 
equation,  see  -for  example  Kanwal  (1971). 

We  need  the  following  lemma. 


Lemma . Let  f(u)  be  a strictly  increasing  and 

d ifferent i ab 1 e function  in  [a , b]  then 


B 


( 1 s _ f f^(u)du 

(a’  “a)  " J[f(u)-f(a)]Q;[f(b)-f(u)]1-a 


0<a<l 


Proof  : Make  the  transformation  y = ^ ■ Thus, 

f (b)-f (a) 


f(1)(u)du 


a. 


[f  u-f  (a)  ] a [f  (b)  -f  (u)]  1-a 


y (1-y)  dy  = B ( 1-a , a) 


0 


since  0<a<l  . When  m is  even,  let  m=2k;  we  have  the 
integral  equation 


oo  , , 

s 1 fh(r)  , j,  K21-3) 

g”(r)  =B(|,i(.k-i))  (y— - > d^- 


Using  the  method  given  in  section  3.4  we  get  the  result 


oo 


-d  U-i  _frx  (2k-3)(2k-5)...3.1fg(y)  2 2 

rdr  I ^ 2k  ^ ' r l /"ob-  n)  I . ,2k  — 2 Cy  r ) dy 


B(i,i(2k-1))  Jy 


that  is 


(\  oo  , 

— d l /-  \ (k-1)  !0j;— i f h(y)  f 2 2\  ~ 2 i 

Fd?J  §2* ( r)  — f — 2 J ^2*32  (y  -r  ) dy 
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Multiply  both  sides  of  this  equality  by  . and  then 

■\r2-z2 

integrate  with  respect  to  r to  get 


oo 

f — 2r 


C(-r4-)*+1g2fc(r')]d] 


z 


oo  oo 


_ (k  2r  h (y)  f.,7.  „2n  3 

* ] n — 2 


,2*-2 


(y2-r2)  2dy  dr. 


Interchanging  the  order  of  integration  on  the  right  hand 
side  of  this  equality  we  obtain 


2r 


r2-z2  LV  rdr 


t(-^)‘+1g!t(r)  ]d. 


oo  y 

_(k-l)!ot_xfh(y) 


if h(y)  f 

iy2*’!  r’-=J)  (y!-r!) 


drdy . 


Since  B ^ , 5 j = tt , we  can  write  the  above  equality 


as 


op  oo 

f T[(-F^)*'1g2t(i')]dr  = (k-1)  !2fc-1[^ 


z(r2-z2)2 


Differentiating  both  sides  of  this  equality  with  respect  to 
z , we  get 


00 

h(z)  = gr-Ci) ! ^Cl  7T7J[(_?37)‘“‘s“(r)]dr] 


z ( r2-z2) 2 


Replacing  2k  by  m we  get 


h(z)  = 


,m  — 2 


?(lm) 


OO 

z z 


2r 


/ 2 2\  2 
(r  -z  ) 
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Even  in  situations  where  h(z)  is  known,  obtaining  h(z) 
by  this  integral  equation  can  be  a difficult  task  for  even 
values  of  m. 


There  are  some  spherical  random  vectors  for  which  the 
probability  density  function  of  the  length  is  known.  We 
can  check  the  result  of  section  4.3  and  section  4.4  by 
reproducing  the  already  known  results  using  the  methods  of 
these  two  sections.  We  begin  with  the  multivariate  normal 
d i st r i but  ion  . 

Let  X be  a random  vector  such  that  X~N(0,  I)  then  by  the 
result  of  the  previous  section  the  probability  density 
function,  say  h(r) , of  R=|X|  satisfies  the  integral  equation 


Assume  that  m is  odd,  for  convenience  let  m=2k+l  , k=l  , 
2, . . . , then  equation  (4.5.1)  can  be  written  as 


4 . 5 The  Multivariate  Normal  Distribution 


(4.5.1) 


Differentiating  both  sides  of  this  equality  with  respect  to 


r and  then  dividing  by  -2r,  we  get 


dy . 


-89- 


Repeat  this  procedure  a further  (k-2)  time  to  get 


1 ^.2 


1_  [2  "5r  _ 2f(k)  fh(y) 

“ B(i,k)  J y2*"1 


dy. 


Differentiating  with  respect  to  r,  we  get 

2r(k)  h (r) 


)fc— 1 NTT 


= re 


-lr2 
2 1 


B(i,k)  r2i_1  ’ 


which  can  be  written  as 


h ( r) =2- 


-ir2 

r2*e  2 


2|(2fc+Dr(21of 1) 


0<r<oo  , 


which  is  the  probability  density  function  of  a X2k+i 
variable . 

We  now  examine  m as  an  even  number.  Assume  that 
k=l,2,...,  then  by  the  discussion  above 


oo 


2 „4r2  _ (2k-3)(2k-5)...l'f  h(y)  i 
' B(J,J(2k-l))  J.J7T72  y2k~2  y 


2k  f (k) 


f h(y) 

1 

dy . 

Ny!-r’  : 

y2k~2 

(4.4)  we  can 

write 

Z)  -J2 

OO 

d f 

t . e 

b — 2 N’T 

dz  J 

z 

dt 


random 


m=2k , 


To  evaluate  the  integral  above  make  the  transformation 
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2 2 2 • • 

s =t  -z  giving  sds  = tdt  and  the  integral  becomes 


-|z 2 °?  -is2 
e 2 e 2 ds 

0 


42t  e 2 


-W2 


we  write  h (z)  as 


h (z)  =2  s!*.1  e 'Z  , 
2*T(k) 


or  the  density  function  of  R 


as 


h (r)  =2  r2t  1 e ^ 


2*T(k) 


0<r<oo , 


which  is  the  probability  density  of  a \ik  random  variable, 


4 . 6 Cauchy  Distribution 

We  can  further  illustrate  the  use  of  the  results  of 
section  4.3  and  section  4.4  by  considering  the  multivariate 
Cauchy  distribution.  The  density  function  of  the  length  of 
the  sum  of  independent  Cauchy  random  vectors  has  been 
obtained  by  other  methods;  thus,  we  can  compare  the  results 
of  this  section  with  other  known  results. 

Let  Z be  a random  vector  and  let  the  probability 
density  function  given  by 


f(z)  = 


r (*»+*) 


(rd)) 


i(m+l) 


(l+z'z) 


\ (m+0 
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then  Z is  said  "to  have  a multivariate  Cauchy  d i str i but i on . 
Let  Zj , Z2 , . . . ,Zn  be  n independent;  Cauchy  random  vectors, 
and  let 

a = [(£  z,)'(t  z,)] 1 

the  distribution  of  R when  n=l  is  known.  We  now  give  a 
derivation  of  the  probability  density  function  of  R when 
n>l . First,  we  state  and  prove  a result  from  complex 
analys i s . 

Jordan’s  Lemma : 

Let  r be  a semicircle,  centered  at  the  origin  and  of 
radius  R,  and  f(z)  be  subject  to  the  conditions 

(i)  f (z)  is  meromorphic  in  the  upper  half  plane. 

(ii)  f(z)->0  uniformly  as  | z | ->oo  for  0<arg  z<tt 

(iii)  m is  positive,  then 

emizf  (z)  dz— tO  as  R— »oo. 

r 

Proof : Phillips  (1956)  page  123. 

By  (ii)  of  the  lemma,  if  R is  sufficiently  large,  we 
have,  for  all  points  on  T |f(z)|<e.  Now 

| exp  miz|  = | exp  miR{cos<?+isin0}  | = exp  -mRsin#. 


Hence , 
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f(z)em<z  dz|  = | [f  (z)emiz  Rei9  id9 


r 0 

7 r 

<e|e_mR  Sln  9 Rd6> 

0 


=2Re 


e-  mR  sin  0d9' 


0 


Now  it;  can  be  proved,  by  considering  "the  sign  of  its 

derivative  or  otherwise,  that  —p9  steadily  decreases  from 

u 

1 to  2/ 7r  as  9 increase  from  0 to  . Hence,  if  O<0<^7r,  then 

sinfl  ^ 2 
9 - * • 

Hence 

7T 

| Jf  (z)em,zdz  | < 2Re  ^e~2mR9 itdO 
r o 

— IA  ( 1 — mR  \ ^7T€ 

— m 1 ^ m • 


from  which  the  lemma  follows. 


This  lemma,  together  with  the  Cauchy  residue,  theorem 
is  useful  in  evaluating  integrals  of  the  form 


, nr  ^ N(z) 

where  Q(z)=p^ 

(i)  D(z)has  no 

that  of  N(z)  by 


and  N(z)  and  D(z)  are  polynomials  such  that 
real  roots;  ( i i ) the  degree  of  D(z)  exceeds 
at  least  one  and  (iii)  m>0 . 
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Consider  "the  "following  integral 


a 
7 r 


-oo 


e,,x 

a2fx2 


dx . 


When  a=l  "this  is  "the  characteristic  function  of 
univariate  Cauchy  density  function.  To  evaluate 
integral,  we  take  the  contour  which  consists  of  the 
segment  [~R,R]  and  a semicircle,  T,  of  radius  R in 
upper  half  of  the  (x,y)  plane. 

Consider  the  following  contour  integral 


e‘,z  . 
dz 


2 , 2 
z +a 


where  the  contour  is  as  described  above.  The  poles  of 
integrand  are  the  roots  of  (z2fa2)=0,  that  is  z=±ai  , 
only  z=ai  lies  within  the  contour.  The  residue  at  z=ai 


isZ 

lim  (z-ai)7 7 

z-»ai  y (z-a l ) (z+ai  ) 


_ e 


— sdb 


2ai 


Thus  the  residue  is  — . 


Observe  that 


I e*8Z  dz  - 
z2+a2  “ 


R 

i 

-R 


Di>x 

e dx  + 


2 , 2 
x fa 


af  e»z 
a2fz2 


dz  = e ' 


the 
th  i s 
1 ine 
the 


the 
but 
i s 


By  Jordan’s  lemma  the  integral  around  T approaches  zero. 
Letting  R tend  to  infinity  we  get 
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oo 

a 
7 r 

-oo 


e”x 

a2+x2 


dx 


e 


— «a 


s>0  . 


The  characteristic  function,  </>(s),  of  the  sum  of  n 
independent  identically  distributed  Cauchy  random  variables 
is  obtained  on  setting  a=n ; giving 
*(s)  = e-ns . 


The  probability  density  function  of  the  sum  of  n 
independent  univariate  Cauchy  random  variates  is  given  by 

f(r)  = * (n2  l r2) 

and  the  probability  density  function  of  the  absolute  value 
of  a single  element  of  the  sum  of  n independent  Cauchy 
random  vectors  is  given  by 


hi(r)  = 2n 


. n2+r2 


The  techniques  developed  in  this  chapter  yield  the 
probability  density  function  of  R= 
number  of  dimensions.  In  the  instance  of  2m+l  dimension, 
we  have  the  integral  equation 


for  an  odd 

»=i  i=i 


oo 


2n 
T n2+r2 


1_  = 2fW!)  (z2~  r2)  m~1 


z2m+i  B(i,m) 


dz  , 


Differentiating  both  sides  with  respect  to  r we  get 
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n -2r 
' (n2+r2)2 


OO 

•I 


2 „2\m-2 


-2r (m- 1 ) r h2m+1(z)  (z  -r  ) 


z 


2m  — 1 


dz 


that:  is 


" (n2+r2) 2 B(l,m) 


OO 


2 J\m  — 2 


(m- 1 ) f h2m+1(z) (z  -r  ) 


2m  — 1 


On  differentiating  a further  (m-2)  times  we  get: 


*■  (n2+r2)m  B ( 


OO 

1 [ 

5’m)  J 


h2m+l(Z) 


,2m  — 1 


dz  . 


D iTTerent i at i ng  with  respect  to  r tor  the  rfinal  time  we  get 


n 


m . 2r 


K ' r 2 , 2\  m+1 

(n  +r  ) 


h2m+l(r) 

r2m_1  B(i  m)’ 


that 


h2m+i(r)  = gSrgT.-2,rr^m+i 


,2m 

/"  2 , 2 \ i 

(n  +r  ) 


_2nm  r(m)  FG)  r2m 


T r(m+|;  '(n2+r2)m+1 


2nr(m+l)  r2m 


r(m+i;  r4/(n2+r2)m+1 


0<r<oo  , 


We  now  state  the  theorem. 


Theorem  14.  Let  X1,X2,...,Xn  be  independent,  identically 
distributed  random  vectors  with  Xx  having  a (2m+l) 
dimensional  Cauchy  distribution. 
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Let 


R2  =]  (ZZ'i)  (EX,) 


then  "the  density  function  of  R is  given  by 


h(r)  = „,,2"  , 


,2m 


6(35 mf 2)  (n2+r2) 


2 , „2\  m+l 


0<r<oo , 


For  the  2m+2  dimensional  case  we  have  the  integral 
equat i on 


00 


n 1 


* (n2+r2) 


h2m+2(z)  (z2-r2)2 


\ (2m  — l) 


z2m  B(i,i(2m+1)) 


dz  , 


D itferent  i at  i ng  “this  expression  m-times  with  respect:  to  r 

we  get: 


00 


n 


(n2+r2)m+1 


h-n*iZ)  (-’-r2)"1  dz 


giving,  by  the  discussion  in  section  4.4, 


^2m+2  ( r ) __  -n  d 


00 


_ f z 1 

'2m  * dr  j (n2+Z2)m+1’>|(z2-r2) 


dz  . 


Make  “the  substitution  that  (z2-r2)=t2,  then 


00 

I? 


/ 2 , 2\  m+1  n 2 2\ 

(n  +z  ) ^(z  -r  ) 


dz 


00 

= 1 
0 


dt 


(n2+t2+r2)m+1 


dt 


Make  the  second  transformation  t=  J(n2+r2)  tan#,  which  gives 
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oo 


dt 


2,^2  , „2\  m+1 


0 


(n2+t2+r2) 


= {n$?Y'm*V  <!0s2me  dl> 
x 7 0 


\ 2m-f~l 

.2-7^2  M B(y(2m+1)) 


n +r 


Thus , we  have 


h2m+2  0)  _ n 


^ = -#(2m+l)  B(y(2m+l))^;f-^p 


= #(2m+l)B(i,i  (2m+l)) 


( 2m-i-3  ) 
2 


giving 


h ( r = 2n r ^ 

2m+2l  7 B(i  m+1)  (2m+3)- 

2 /-  2 , 2\  2 

(n  +r  ) 

The  probability  density  function  is  thus  given  for  all 
value  of  n>l . 

We  could  have  obtained  these  results  by  noting  that  an 
m-d imens ional  spherical  Cauchy  random  vector  could  be 
obtained  by  considering  Zx/Zm+1 , Z2/Zm+1 , . . . , Zm/Zm+1  where  Z, 
i=l , 2 , . . . , m+1  are  independent  N(0,1)  random  variables. 
Then  ( Z1/Zm_|_1 , Zj/Z^^j  , . . . ,Zm/Zm+1)  has  an  m-variate  Cauchy 

distribution.  Earlier,  we  showed  by  the  use  of 

characteristic  function  that  the  sum  of  independent  Cauchy 
random  variables  has  a Cauchy  distribution. ; If  we,  in 
fact  denote  this  random  variable  by  X then  the  probability 
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density  function  of  Y=X/n  would  be  given  by 


f(y)=  i 


It  follows  that  the  length  of  the  sum  of  independent  m- 
variate  Cauchy  random  vector  has  the  same  distribution  as 
Xm/riXx  where  xl  and  Xm  are  independent  chi  square  random 
variables  with  1 and  m degrees  of  freedom  respectively. 


4 . 7 Uniform  Distribution  on  the  Unit  Sphere 

In  Chapter  III,  we  gave  a detailed  discussion  of  the 
sum  of  independent  random  vectors,  each  vector  having  a 
uniform  distribution  on  the  unit  sphere.  There  we 
discussed  the  case  of  an  odd  number  of  dimensions.  In 
Chapter  IV,  we  gave  a method  which  allows  us  to  produce 
results  whether  m is  even  or  odd. 

When  m is  even  we  have  that  ^ ’ 2 ’ ' • • Thus,  for 

independent,  identically  distributed  random  variables  W-l  , 
W2 , . . . , W n with  W1~W^,  we  must  determine  the  probability 

n 

density  of  J^W,.  This  is  not  an  easy  task  when  m is  even, 

i=i 

since  the  density  of  is  no  longer  a polynomial  of  finite 
degree.  This  prevents  us  from  obtaining  the  density 

n 

function  of  ^W,-  using  the  convolution  theorem.  Thus,  the 

i—  1 

method  in  section  4.4  is  useful  only  with  odd  values  of  m. 
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4 . 8 The  Uniform  D i str  i but:  i on  on  the  M-bal 1 


We  can  develop  some  results  for  a random  vector  having 

a uniform  distribution  on  the  unit  m-bal 1 by  using  some  of 

the  results  which  we  proved  for  a random  vector  having  a 

uniform  distribution  on  the  unit  m-sphere . We  can  say  that 

a random  vector  B has  a uniform  distribution  on  the  unit  111- 

ball  and  write  B~U(Bm)  if  P(||B||<r)  = = rm  whenever 

Vm  (.11 

D 

0<r<l  and  jf=-r.~U (Sm)  , where  Vm(t)  denotes  the  volume  of  a m- 

lel 

dimensional  ball  of  radius  t. 


Theorem  15 : Let  U be  a random  vector  having  a uniform 

distribution  on  the  unit  m-sphere  and  let  A be  a random 
variable  such  that  A~B(1 ,m)  and  independent  of  U,  then 

B = A U~U (Bm) . 


Proof : The  transformation  Y=aU  maps  a point  on  the  surface 
of  the  unit  m-sphere  to  a point  on  the  surface  of  the  m- 
sphere  of  radius  a.  The  surface  area  of  the  m-sphere  of 
radius  a is 


2a™-1(r(j))™ 

r(!m) 


Thus,  the  conditional  density  of  AU  is  given  by 


f (a  U | A — a) 


r(» 

2am_1  (f(|))  m 
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We  obtain  the  density  of  B = AU  by  multiplying  conditional 
density  given  above  by  the  marginal  density  of  A and  then 
integrating  the  expression  with  respect  to  a.  Thus 


g(B) 


" f (§m)  m am-1 
l 2am_1  (f(j)  )m 


1 

0 


mT  (lm) 

2(r  (hr 


da 


mf  (|m) 

2(r(i))n 


Since  0<A<1 , multiplication  of  U by  A maps  onto  the  unit 
ball  giving 


f(B) 


mr(|m) 

2(r(i))" 


Theorem  16:  Let  B be  a random  vector  having  a uniform 
distribution  on  the  unit  ball,  and  let  a1  be  a unit  vector 
then 

a B,~W"+1. 

Proof : There  are  several  methods  which  we  can  use  to  prove 
the  result.  However,  our  proof  is  based  on  the  results 
which  we  used  in  our  discussion.  Let 


<£(t)  = g elta^ 
iZ _ 

m 


i tAa/.r^Ti 


=8  e 


, (where  Z ~N (0 , !))• 
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Thus 


<f>( t)  =8  e11^  = Jmam_1Vj,  (at)  da 


=|mam_1r(j/+l)  (|at)  ~1'  Jj,(at)  d 


a 


= (|)_I/2r(i/+2)  a^^^Cat)  da. 

0 

On  page  373,  Watson  (1944)  states  that 


Ju+J,+i(z)  - 


v+i 


2I/r(^+i). 


Ju  (rs  i n0)  S i nu+10  cos2I/+10  d 9 


0 


which  is  valid  when  the  real  part  oT  v exceeds  -1. 

Make  the  transformation  t=sin0  and  then  put  i/=0  to  get 


J. 

Jau+1  Ju  (at ) da  = t-1Ju+1(t). 


Thus  we  have 


*(t)=S  eitA^ 


-(It)-12  2T (r+2)  t J^+1(t) 


:(It)_(l/+l)r(r+2)JJ/+1(t) 


=V|/+i(t)  =8  e 


itW 


v+i 


where  W~W 
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In  other  words,  a single  element  of  a random  vector 
having  a uniform  distribution  on  the  unit  m-bal 1 has  the 
same  distribution  as  a single  element  of  a random  vector 
having  a uniform  distribution  on  the  unit  (m+2) -sphere . 

We  have  established  the  relationship  between  h2m^_1(r) 
and  h^r)  for  a spherical  random  vector,  where  h2m^_1(r)  is 
the  probability  density  function  of  the  resultant  of  the 
(2m+l) -component  vector  and  h^r)  is  the  probability 
density  function  of  the  absolute  value  of  a single  element 
of  that  vector.  For  the  (2m+l)  dimensional  unit  sphere 
hj(r)  is  the  density  of 


Y 


n 


iEw«i 


9 


f(2m+l)-l 

where  W^~W  and  W!,W2,  • • . ,Wn  are  independent;  whereas 

for  the  (2m+l)  dimensional  unit  ball 


Y 


iE 

V — 1 


W,l 


i(2m+3)-l 


where  W;~W2 


The  results 
applicable  to  the 


of  section  1.4  and  section  1.2 
case  of  independent  random  vectors, 


are 

each 


of  which  has  a uniform  distribution  on  the  ball. 


CHAPTER  V . 

SUMMARY  AND  FUTURE  RESEARCH 
5 . 1 Summary 


In  Chapter  I we  gave  some  of  the  properties  of  a random 
vector,  which  is  uniformly  distributed  on  the  unit  m- 
sphere . We  discussed  some  results  which  were  given  by  Saw 
on  the  distribution  of  Rn , the  length  of  the  sum  of  n 
independent  vectors,  each  of  which  has  a uniform 
distribution  on  the  unit  m-sphere . 

In  Chapter  II  we  dealt  with  the  computational  problems 
encountered  in  using  the  procedure  given  by  Saw.  A part  of 
our  investigation  was  to  approximate  the  distribution  of  Rn 
using  the  probability  density  functions  of  a set  of  Beta 
random  variables.  Our  results  though  encouraging  did  not 
improve  on  the  method  used  by  Saw. 

In  Chapter  III  we  gave  a result  for  the  probability 
density  of  the  length  of  the  sum  of  independent, 
identically  distributed  random  vectors,  each  of  which  have 
a uniform  distribution  on  the  unit  m-sphere;  here  m>3  is 
an  odd  integer.  In  Chapter  IV  we  presented  generalizations 
of  the  earlier  chapters.  We  considered  spherical 
random  vector  and  derived  probability  density  functions  of 
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their  length  in  terms  of  the  probability  density  functions 
of  the  absolute  value  of  a single  element.  We  also 
described  a method  for  finding  the  density  of  the  length 
for  an  odd  number  of  dimensions  and  we  were  able  in  some 
situations  to  obtain  results  for  even  number  dimensions. 

5 . 2 Future  Research . 

In  Chapter  II  a failure  in  a part  of  this  study  was 
reported.  Further  investigation  of  the  problem  is 
necessary,  and  it  is  my  intention  to  return  to  this 
problem . 

Another  area  of  interest  is  that  of  elliptical  random 
variables.  In  the  same  way  that  the  family  of  spherical 
random  variables  contains  X~N(0,I)  as  a special  case,  the 
family  of  elliptical  random  variables  contains  X~N(0,E)  as 
a special  case.  Many  of  the  results  which  have  been  proven 
for  a multivariate  normal  distribution  are  also  true  for  a 
spherical  random  vector.  Although  the  mathematics  required 
for  an  elliptical  random  vector  are  somewhat  more 
difficult,  it  is  my  belief  that  we  can  establish  some 
results  for  elliptical  random  vectors  which  are  similar  to 
results  for  X~N(0,E).  It  is  my  intention  to  pursue 
research  along  these  lines. 
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